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PRESENTACION 


Este ejemplar dela Revista SELPER forma parte de 
dd a casta icararos ca han co lr: 
“gran esfuerzo, en un Plan Editorial de 
Emergencia, destinado a recupera la continuidad 
editorial de nuestra Sociedad, iniciacia durante la 
Sede de SELPER en Chile (1983-86), la que, lamenta- 
blemente e 06 se pudo 

juir en la forma en que hubiésemos 

E E leones che 
dos simbólicamente como Septiembre 1986, Sep- 
tiembre 1987, Septiembre 1988 y Septiembre 1989, 
para representar la continuidad en esos años y co- 
menzar a publicar normalmente en 1990. Dado. 
son Revistas editadas con posterioridad (St "e 
1989, fecha de materialización del Pian Editorial de 
Emergencia), no e han incluido Secciones de Not 
cias nÍ Actividades Iberoamericanas e Intemnaciona- 
les, las que sí aparecerán en los números regulares 
qué serán publicados durante 1890 en forma semes- 
tral y, si ol iberoamericano e internacional es 
adecuado, en forma trimestral. 

Los artículos técnicos seleccionados para la Revis- 
ta SELPER 1996, corresy atrabajos presenta- 
dos en el Seminario de ja Regionalcoor 
zado porla Unión Intemacional de 
¿as (US), UNESCO, SÉLPER y el Inettuto de Pas" 

quis Engacas (INE) de Brasil (Sao José dos 

dampos, 8P, Bras, Diciembre 1989), La cortinas, 
ción de estos artículos técnicos la pane 
medular de la Revista SELPER 1987, Para la Revista 
SELPER 1988 se ha seleccionado 
completo artículo del Dr. 
¡Nacional Autónoma de México, 
torización Matemática cie Imágenes 











sena creido conveniente 

mido jemplficador trabajo del Ing. 

Camy lic sión Nacional de Investigaciones Espa 
ciales de Argertina, CNIE), destinado a profesionales 





no tan especializados, Finalmente, para la Revista 
SELPER 1989, se continuó con el siguiente capítulo 
asco por el Dr. Jorgo Lira (UNAM) ALsotr Table 

spera: 





maciones Reversibles y Realces 
mos que estos temas sean de ner de yrecibira fro, 
para las próximas Revistas SELPER, otras valiosas 
contribuciones de la comunidad profesional latinoa:- 
moricana y mundial. 
Probablemente lame la atención alpor qué de eso 
aa se materializa, 
asrazones son constituyen un compro: 
mio de Poner, de comme y coplas periodo. 
amino de las acividades de SELPER a lola 
larotación de sus Sedes. Durante laSede de SELPER 
Brasi(1986-89) se obtuvieron inmensos logrosen 
amplar a cooperación intemacional, organización 
imposios Latinoamericanos en diversos países 
de Nudata Amibca en cooraracón ón la Sedo 
SELPER en Brasil, reconocimientos intemacionales 
hacía nuestra Sociedad y diversas actividades de 
apoyo dento y fuera de Iberoamérica. Todo ello 
demuestra la madurez y consolidación que ha logra- 
GOSELPER, Lamentablemente, por razones de uer. 
a A 
Conscientes dela importancia que tienen las Publ 
caciones como agentes de comunicación y pertec- 
cionamiento permanente de nuestros Socios, ya que 
no es posible para todos en los diversos 
Cursos de Perfeccionamiento, os, Simpo- 











sios y otras actividades, tadas O co-patrocina. 
das por SELPER y, conla misma visión y generosidad 
¡con que la Sede SELPER Brasil abordó la realización 
de Reuniones Plenarias y Simposios Latinoamerica- 
"nos fuera desu país, también haexistido una apertura 
Para realza un Plan Estoril a más largo plazo, 
raspasando importantes responsabil 
dades al Comas Exitodal que tantos éxhos logr 
reposo O aj 
yde Arbitraje a 
vel , Compuestos por profesionales 
de probada ascelencay scicaciónaSELPE, o 
ala esperada decisión de las futuras Sedes de conti 
inuar con un Comité Editorial más estable indepen- 
¿ientemente de la rotación de éstas. deberá permitir 
lograr la tan ansiaca continuidad editorial que tanto 
necesita nuestra Sociedad para fortalecer y muliplk 
Car sún más los áxics cosechados en todos estos 
años, desde sus inicios en Ecuador, pasando 
ia a Argenta y pomos nos pss 


lla valiosa colaboración presta- 

O os 

Nacionales SELPER, en icular los de 

Brasil, Chile, Colombia, Ecuador, México, 

Ce e hermanos de ql Y 

Ad E 

E ca nenas 

Igualeecie alos Socios Inuciones yEmeresns de 

Honorarios de Canadá, 


los Capítulos 
España, EE JU. Francs, Holanda y Portal [gua 
mente, ls organizaciones Ineracionales como 


'uropea (ESA), Naciones Unid 
PoxiON (SO USS COSPANAGSSy dia ¿que tan dect ca 


también han aportado a los l0 SEL- 
PER, Pialmento, pata la realzación do esto Plan 
Estoral de Emergencia, vayan, nuestros encaros 
reconocimientos Publicaciones SEL- 
PER Chile, Ecitorial Universitaria, Impresora América 
y, Por supuesto, al Instituto Geográfico Mitar (GM) 

1 a su Director, Brigadier Don 
ovario Casilo Vel subdirector Téshico TEL Lada. 
ro Rivas, TCL Sergio Matus y MY. Mario Molina, pue 
esa colaboración será vital para consolidar est 
Editorial a largo plazo de SELPER. 























para SELPER Chile, un 
¿ran honor poder comimuar cortbuyendo on Forma 
tan decisiva a la mayor consolidación de SELPER; y, 
para SELPER Argentina, una gran satisfacción y un 
desafo el poder continuar avanzando cada vez más 
en estas actividades en beneficio de SELPER y sus 
ideales de cooperación y hermandad iberoamer 

y universal. 


Sedo SELPER Bras (1986-1989) 
PROF. MAURICIO ARAYA-FIQUEROA 
Vice Presidente (Ecuador, 1980-89) 
Presione (Cra, 1989.26) 
Vice-Presidene y Drecior Era! 
¡Sede SELPER Brasa (19069) 

ng. MIGUEL SANCHEZ-PEÑA. 

Vice Presente SELPER 

Sedo SELPER Bras (1986-89) 
residente Comité Organizador. 

Sede SELPER Argentina (1989.91) 


EDITORIAL 


Thisissue, thatis par ofa series offour (4) SELPER 
Reviews, constíttes a great effort to continue with the 
Editorial Programme that was begun by SELPER- 
Chile during 1983-26. Unfortunately and due to ser- 
0us problems, this could not be continued. Because 
otthis, we are publishing now these editions 
which wil represent the continulty along years 1986- 
1989. These issues have been symbolically dated 
September 1986, September 1987, September 1988 
September 1989. Due to the fact nat these Re- 
views were edited in September 1989 (date of the 
materialization of this Emergency Editorial Plan), nai- 
ther international and Iberoamerican Activities nor 
Nows were included, These Activities and News aro 
going to appear in the normal sections in the regular 
issues which wi bo during 1990, every six 
months. but íf the Iberomerican and International 
duppe/ tsutabl, SELPER Review wl published 
quanert 





The technical papers selected for SELPER Review 
1986, come from the Seminar on Regional Geology, 
coordanize by interatonal Union ol Geological 

1ces (IUGS), UNESCO, SELPER and INPE 
(Brasilan Space Research Institute), heldin Sao Jose 
dos Campos, Brazll, on December 1985, SELPER 
Review 1987 was shaped with most of these papers 

ible at the Editorial Committeo. An important 
paper about Mathematical Charactorization of Digital 
Images by Dr. Jorge Lia (UNAM, NatignalU 
lindoo) maños SELPEN Rewiow 1988 ogainar uh 
a surmarzes paper by Eng, Marcelo Cami CNE, 
Argentinian Spaco Research Commission). Finally, 
SÉLPER Review 1989 contains the next chapter wr: 
100 by Dr. Jorge Lira (UNAM) about Reversible Trans- 
formations and Global Enhancementes. We hope 
those papers wil be of interest and hope to recere 
'morecontributionsforthenext SELPER Reviews from 
Latin American and worldwide community in the nexd 
futuro, 

The reasons to do this Emergency Editorial Pian, at 
this time, constitute a commiiment and a constant 
improvement n the activities of SELPER all along the 
rotation of ts Seats. In Brazil (during 1986-89), 
SELPER got alot of achievements in the international 
Opening and coorganization of Latin-American Sym- 
posia «which took place in different Latin American 
countries. Our Society also received addticnal iner- 
ational recognition, as well as international suppon. 
The above showed that our Society had reached an 
important consolidation and maturiy. In spite ofthist 
was, Unfortunately, not possible for our Society to 
reinforce SELPER's Editing activities, 

Aware of the importance that these publications 
haveas agents ofcommunicationand constant tooo! 
improvementfor our members, sinceitis not possible 
for all of them to take part in the diferent training 
courses, seminars, symposia or anyother activities 
which are organized or co-sponsored by SELPER 
and, wih the same view and ty which 
SELPER/Brazil organized Plenary Meetings and 
Latin-American Symposia with abroad there has also 
xisted the same will to accept and to create a long- 
term Edition Program by applying and transtering 




















Sepiente 1996 SELPER 


important and significant responsibiitias to the Edit- 
ing Commíttes that had so successful participation at 
he SELPER Seat in Chile (1983-86). lts, plus the 
integration of both a Specialists Commtteo and Roto: 
reo! thatare consttuted by professionals of 
high excellence and devotion to SELPER, together 
wth the expected decision of the futuro SELPER. 
Seatsin orderto continue working with a more perma: 
nentEdting Committes thatis not goingtobe affected 
by SELPER Sears rotation, we hope that it wil alowto 
getthe so expected editing continulty that our society 


needs to reinforce and even more tho 
Achievements got in the past, since its begining in 
Ecuador, going through Chile, Brazil, Argentina and 


Other brother countries in the near futuro. 
We also want to thank the valuable cooperation 
given by the diflerent members of SELPER and the 
Institutions where the National Chapters of SELPER 
worked, specially hose from Argentina, Brazil, Chio, 
Colombia, Ecuador, México, Perú and other brother 
countries of South America and Central America, 
whose suppor and cooperation and permanent par: 
'icipation were fundamental forthe creation and exis- 
once of SELPER, We also want to thank to all tho 
Members, Enterprises and Institutions ofthe SELPER 
Honorary Chapters, specialy those of Canada, 
France, . Holland, Portugal, Spain and 
USA. Also to those international organizations like 
European Space Agency (ESA), United Nations, 
1UGS, COSPAR/ICSU and others who have been so 
decisivo for SELPER success. Last but no loast, wo 
want to thank the valuable support given by tha 
'SELPER-Chile Editing Commites, Editoral Universi 
tara, Impresora América and Instituto Geográfico 
Miltar (1GM) from Chile - specially to its Director, 
Brigadier General Mr. Roberto Castilo V,; Technical 
Subdirector, LTC Lautaro Rivas; LTC Sergio Matus 
¡and MY, Mario Molinawhose colaboration will bo of 
real importance to pertorm this SELPER's long:term 
Edition 
ln this way, ít is a real proud for SELPER Brazil to 
offer an updated ecition programm at the end of ts 
mandate; ts a real proud for SELPER-Chile to con- 
tine! in such a decisive way lo the greater 
Consolidation of SELPER; ít is a real challenge and 
satistactión for SELPER- to continue im- 
these activities for the benefit oPSELPER and 
important aims of Iberoamerican and Universal 
brotherhood and cooperation. 








Prot. MAURICIO ARAYA-FIGUEROA 
'Vice-President (Ecuador, 1980-83) 
President (Chile, 1983-86) 

en preident and ESieno Director 
SELPER Seat Brazil (1986-29) 


Ing. MIGUEL SANCHEZ-PEÑA 
Ver eresdent ol SELPEn 
Beso! Orar Cominitos 
SELPER Seat Argentina (1989-91) 


SELPER seentre 1988 


ANTECEDENTES 
DE SELPER 


LaSociedadde Especialistas Latinoamericanos en 
Percepción Remota (SELPER) es una organización 
Técnico-Cientifica sin fines de lucro, Los principales 
objetivos de SELPER son promover el perfecciona: 
miento de sus asociados y apoyar la 
eficaz entre diferentes países, con la finalidad de 
¡provecharal máximola capacidad existente eincen- 
isa a integración entre los mismos a través dela 
Cienciay la Tecnología. Estos objetivos permiten una 
integración a nivel nacional, regional continental e 
interregional (Intercontinental) uniendo. naciones 
con problemas similares. 








'SELPER fue creada en Quito, Ecuador, en noviem- 
bre de 1980, durante la Primera Reunión de Especia- 
Iistas en Percepción Remota. En este sentido fue 
cumplida una aspiración de latinoame- 
fcanos,los cualesno contaban conuna organización 
para ol desarrolo de programas de integración 
cooperación en América Latina. Esta ir ración, O. 

inalmente a nivel individual, fue extendida a nivel 

tucional y hoy varias organizaciones intemacio- 
nales participan directa o indirectamente de las activ. 
dades de SELPER. 


Para minimizar costos y para facilitar decisiones, la 
Era ad ca [arc 68 Muy Sim» 
9: (a) Un Consejo Central forme 
[osdolpale designado comosade de SELPEN pero” 
do normal de dos, años) el cual toma las decisiones 
cutivas, lavando en cuenta el mayor beneficio de 
la sociedad y sus ideas. (0) Coordinadores naciona- 
les, formados por profesionales de cada país y repre- 
sentada por un Coordinador Nacional que coordina 
las actividades de SELPER a nivel de su país y 
transmite las aspiraciones de sus miembros a la Sede 
Central de SELPER y/o durante las Reuniones Piena- 
rias de SELPER. (c). iblea General, es organizada 
durante las reuniones plenarias de SEI y está 
compuesta por socios o sus representantes (Coord 
nadores Nacionales). 























La! 





ón aSELPER incluye tres categorías prín- 
SocioIndividual, normalmente profesiona: 
'eon actividad en percepción remota; b) Socio 





le 
Institucional, normalmente organizaciones sin fines 





deluero como Universidades u otras; c) Socios Patro- 
cinadores, normalmente organizaciones comerciales. 
que deseen cooperar con SELPER. 


La Sede de SELPER es rotativa y se ha ubicado en 
Ecuador (1980-83), Chile (1983-86), Brasil (13 
es Argentina (1989-91) y Perú (1991- 

'ELPER realiza o co-patrocina Cursos, Seminarios, 
Simposios y otros eventos, destacando las Reunio- 
nes Plenarias asociadas con Simposios de carácter 
«continental, las que se han desarrollado en Ecuador 

(1980 y 1982), Chile ar 1984), Canadá (1985), 

asil(1986), Colombia(1 Mic (1968). Argen 
tina (1989) y al cumplirse diez (10) años de SELPER, 
Ecuador (1990). La Sociedad realiza diversas publica- 
ciones, como Memorias Anuales (correspondientesa 
'cada Reunión Plenaria y Simposio asociado), Calen- 
darios Técnicos o Almanaques (anual o bianual), 
Revista Técnica (normalmente semestral otimestra) 
y otras publicaciones especiales ocasionales. 








SELPER 
ANTECEDENTS 


he Society of Latin American Specialists on Re- 
mote Sensing (SELPER) is a technical scientific and 
non remunerativo Society, Which is aimed at the 
achievement of the steadlast professional improve- 
'mentof!ts members and promoting.or suppoingthe 
effective cooperation among diferent countries, in 
order to take advantage of the existing capacit 
the extreme limit andlto integrate the countries one to 
each other by means of science and technology. 
These purposes are looking for an integration to Na: 
tonal, Regional, Continental and inter-regional level, 
in this way nations having similar proble 
SELPER was created ín Quito, Ecuador, on n 
ber 1980, during the first meeting of Latin American 
Specialists on Remote Sensing. In this sense and old 
dream was performed by professionals of the region 
bo dd nck have an organization 10 develop sl 
cient, programs of integration and cooperation in 
Iberoamerica. This integration, atthe begining, was to 
an individual level, being afterly 10 Institu- 
tonal level by which several International Organ 
tons are participating ofthe SELPER Activities directly 
or inderecti. 


In ordertoreduce the administrations costs, and to 
enable the decisions, the SELPER administrative 
organization is very simple: a) Main board which 
formed by professionals belonging to the SELPER 
Seat Country. (Normal tíme-period o! two years). The 
Directory takes the executive decisions taking 
account the best benefit of the Society and ita ide: 
b) National Chapters formed by| jonals in every 
havo a representativo that is tho Nation: 

'which coordinates the SELPER activi 
tonationallevelbbesidestransmting the aspirations of 
he members (Nationals) 1o the SELPER Central Seat 
¡or during the SELPER Meeting.) Gene 
Assembly whichisconstituted or organized duringthe 
'SELPER Plenary Meetings, s by Members ofSELPER 
'orts representativos (National Coorcinators) 

Thero are three main diferent categories for 

SELPER affilation: a) Individual Member, usually 

having activities on remote sensing; b) 
Institutional Member, usually entities having non 
remunerativo purpose like Universities and others; c) 
Sponsor Member, usually companies or other institu: 
tons who wish to cooperate with SELPER, 

SELPER Central Seat normaly rotaes every wo 
years: Ecuador (1980-83), Chile (1883-86), Brazl 
(1986-89) and Argentina (1989.91) and Peru (1991- 
$3) in the next future. SELPER organizes or co-spon- 
sor Courses, Seminars, Symposia and other events, 
being the most relevants the Plenary Meeting and 
associated Symposia: Ecuador (1980 and 1982), 
Chile (1983 and 1984), Canada (1985), Brasil (1986), 
Colombia (1987), Mexico (1988). Argentina (1989) 
and Ecuador (1980) to celebrate the Terth Aniversary 
ofSELPER, The Society publishes Annual Proceed- 
ings (telated to Plenary Meetings and associated 
Symposia), Technical Calendar, (annually or every 
¡wo years), Technical Review (usually every six 
months or quarterly) and other special documents or 
books. 
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ACTIVIDADES LATINOAMERICANAS 
ACTIVIDADES INTERNACIONALES 


Estas secciones serán incluídas en los números regulares de la Revista SELPER, a publicarse en 1990. 
Por ahora, conozcamos algo del Centro de Levantamientos Integrados de Recursos Naturales por 
Sensores Remotos (CLIRSEN) del Ecuador, Primera Sede Histórica de SELPER (1980-1983). Actualmente 
CLIRSEN está ampliando su capacidad de recepción y procesamiento para los satéltes Landsat (MSS y 
TM), SPOT, ERS-1 y satóltes meteorológicos. 













CENTRO DE LEVANTAMIENTOS INTEGRADOS DE 
RECURSOS NATURALES POR SENSORES REMOTOS 





PAZ Y MIÑO S/N, EDIFICIO 16M, APARTADO 8216. 
QUITO, ECUADOR. TEL: 540-585. TLX: 2775 CLIRSEN» 


Ampliación de la Estación del 
Cotopaxi para Recepción y 
Procesamiento de Datos 





. 
PRIMERA SEDE SELPER LANDSAT TM 





TIROS 


AREA DE COBERTURA 
DE LA 
ESTACION DEL COTOPAXI 


CLIRSEN e IGM del Ecuador: Sede de la Primera Reunión que 
dio nacimiento a SELPER, Quito, Ecuador, 24-27 Noviembre 1980. 
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CARACTERIZACION MATEMATICA 


DE IMAGENES DIGITALES 





Dr. Jorge Lira 


Unidad de Percepción Remota. Instituto de Geofísica. 
Universidad Nacional Autónoma de México (UNAM). 
ALP. 22-560. Deleg, Tlalpan. 14000 México D.F. MEXICO 





Resumen 


Se entrega una completa revisión sobre este tema, que constituye el 
quinto Capítulo del libro “Introducción al Procesamiento Digital de 
Imágenes”, del mismo autor. Para un mejor tratamiento del tema, éste 
se ha subdividido en: 5.0. Prólogo; 5.1. Consideraciones Generales; 
5.2. Proceso de Medida sobre un CIV; 5.3. Teorema de Muestreo; 5.4. 
Operadores de Caplura de una Escena; 5.5. Caracterización 
Estocástica de una Imagen Multiespectral; 5.5. Referencias. Las Tablas 
y Figuras se han intercalado en el texto, el que ha mantenido su 
numeración original, pues es laintención publicar otros Capítulos sobre 
eltemade modo que sea posible, através dela Revista SELPER, contar 
conla colección completa (o al menos delos temas más específicos) de 
este libro, lo que ocurrirá de acuerdo a la receptividad que esta idea 
pueda tener entre los lectores de la Revista SELPER. 








PROLOGO 


as sondas espaciales enviadas a 
capturar datos de los cuerpos que 
forman nuestro sistema solar mar- 
can una pauta bien definida en el 








desarrollo del Procesamiento Digital 


de Imágenes. La enorme cantidad de medidas 
hechas a distancias relativamente cercanas a la 
superficie de planetas y satáites mostraron el 


potencialdelacolección de datos en formado imáge- 
nes. Debido a esto y hace ya más de quince años, se 
decidió que ctas bien podrían ser procesadas digl 
talmente por computadora, Ahora existen más de 
150.000 imágenes digitales del sistema solar que han 
revelado detalles insospechados y sorprendentes de 
la evolución del conjunto sol-planetas-satélites-aste- 
roides. En la actualidad no es posible concebir las 
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ciencias planetarias sin el procesamiento digital de 
imágenes, como tampoco es factible pensar en la 
meteorología o la Percepción Remota, aplicadas ala 
prospección geofísica, sin el concurso del procesa: 
miento de imágenes. 


El análisis de imágenes se emplea en muchas 
áreas del conocimiento científico y su uso es en 
verdad creciente, En física, las imágenes del micros- 
copio electrónico son estudiadas para obtener las 
propiedades de la superficie de un sólido. La morto- 
logía de fagocitos es determinada en Biomedicina al 
emplear imágenes del microscopio óptico. Las fsu- 
ras, defectos e inhomogeneidades de una cierta pie- 
zamanufacturada enuna/ábrica, pueden seraprecia- 
das manejando convenientemente la neutrografía 
correspondiente adicha pieza. Ahora en Medicina, es 
posible realzar en color una radiografía digital con el 
Objeto de apreciar con mayor facilidad tumores y 
¿tros problemas en el paciente respectivo. Muchas 
aplicaciones más del procesamiento de imágenes 
existen en diversas áreas de la investigación experí- 
imentaltalescomo: Ingenieríade Suelos, Arqueología 
y Astrofísica. 


No puede verse una imagen digital como simple- 
mente una colección de números arregiados en for- 
ma matricial; esto sería empobrecer su enorme po- 
tencial, Una imagenno.es nada más una colección bi- 
imensional de datos, sino más bien la distribución 
espacial de la respuesta espectral de los elementos 
que componen la escena, entencióndose por res- 
puesta espectral aa variación de intensidad de ener- 
¡glaradiante en función de la longitud de onda de esta 
radiación. Recordemos que aún la radiación corpus- 
cular tiene una longitud de onda asociada. Las gráf- 
cassonunasubclase de imágenes y pueden conside- 
rarse también como una distribución espacial del 
valor de una cierta propiedad de un sistema que 
puede o no tener un significado fsico inmediato. Las 
imágenes tienen una propiedad de conjunto que se 
refierealainterrelación que existe entre los elementos 
de la escena, siendo esta interrelación dependiente 
del sistema físico u Objeto matemático que estemos 
estudiando. En otras palabras, una imagen es una 
colección de objetos y regiones aunada a la interrela- 
ción que existe entre ellos, Esto implica un gran 
potencial de información que demanda de avanza- 
das técnicas matemáticas y computacionales para 
Tacita el camino del modelaje referido al compora- 
miento del sistema bajo estudio. 





Para analizar una imagen se requiere, como se 
menciona en el párrafo anterior, del más avanzado. 
conocimiento científico en fisica, matemáticas y com- 
putación. Enla presente obrase ha dividido el análisis 
de una imagen en dos grandes ramas: El Procesa- 
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miento de Imágenes y el Reconocimiento de Patro- 
es de las mismas; aquí se estudia únicamente la 
Primera rama. La primera de éstas, se entiende aquí 
como un conjunto bien definido y clasificado de 
realces, que dejan ista la imagen para una descrp- 
ción matemática de las regiones que la componen. 
Las regiones están formadas a su vez por patrones 
espaciales y espectrales, de tal forma que al recono- 
cer estos patrones por medio de modelos matemáti- 
os, implica la descripción de las regiones que for- 
iman una imagen. En otras palabras, reconocer un 
patrónrequiere primero de su manifestación evidente 
por medio de un procesamiento bien definido a la 
imagen y, segundo, de una descripción cuantitativa 
de sus propiedades morfológicas y estadísticas. 


Enlos primeros cuatro capítulos, que constituyen 
'unanarración y un análisis delos aspectos generales 
de una imagen digital, se introduce al lector a los 
elementos y estructuras que componen una imagen. 

Para aquel que desee estudiar estos conceptos, 
puede consultar el libro: La Percepción Remota, del 
mismo autor. Debido a que existe este libro, no se ha 
considerado necesario publicar primero los mencio- 
nados capitulos, por lo que la presente publicación 
constituye la muestra inicial del material que compo- 
ne el libro: “Introducción al Procesamiento Digital de 
Imágenes”, que elabora actualmente el autor; los 
otros capítulos serán publicados subsecuentemente 
con la intención de recibir los comentarios y su- 
¡gerencias dela comunidad y lograr asíun manuscrito 
final de ata calidad. Por el momento no se incluyen 
imágenes, pues esto elevaría enormemente los cos- 
tosdeedición; enlaversiónfinal se incluirán ejercicios 
y fotografías ustrando los diversos conceptos y 
procesos quese discuten aqui. Las referencias se han 
incluido en el trabajo por medio de números de 
formato pequeño a manera de supraíndices altexto. 


El material que se presenta aquíha sido desarrolla- 
do durante mis años de impartir el curso: Procesa- 
miento Digtal de Imágenes y Reconocimiento de 
Patrones, de la Maestría de Computación del IIMAS- 
'CCH. Elapoyo que he recibido de lacoordinación de 
esta maestría ha sido constante y mis alumnos me 
han ayudado a depurar parte de lo que se discute 
aquí estoy agradecido por esto. Se hace énfasis en 
aplicaciones del procesamiento de imágenes de sa- 
16lte (Percepción Remota) pues ésta ha sido el área 
principal de investigación del autor. Sin embargo, los 
conceptos y métodos desarrollados son de alcance 
general. El procesamiento digital de imágenes es 
cada vez más empleado en universidades, institutos 
de investigación, industrias, hosptales y tecnológ- 
eos, porlo que espero y deseo que esta obra sea de 
utiidad e interés para la comunidad dedicada al 
análisis de imágenes. 
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5, 1.- CARACTERIZACION MATEMATICA 
DE IMAGENES 
MULTIESPECTRALES 


a cuantificación de la información multiespectral 
establece las bases para la modelación del sistema 
físico bajo estudio y deja el camino abierto a nuevas 
técnicas de análisis de imágenes cigitales por medio 
de algoritmos computacionales. 


5.1.1. Consideraciones Generales 





Enfrente de una pantalla de despliegue blanco y 
negro, se encuentra un radiólogo estudiando una 
radiografía de un posible tumor cancerígeno, con 
base en la interpretación que él haga de esta informa- 
ción, se decidirá si se procede a una intervención 
Quirúrgica o no. La decisión sin embargo se presenta 
difícilpues laradiografa tiene poco contraste yl aka 
definición, Un físico está analizando las primeras 
fotografías de un plasma generado por un reactor 
nuclear de fusión, el científico cree observar las inos- 
tabilidades predichas por la teoría de la fusión nu- 
clear; desafortunadamente, estos datos están conta- 
minados por el ruido inherente al arreglo experimen- 
tal correspondiente y por lo delicado de las medidas 
hechas, de aquí que se complique la interpretación 
Que el físico desea hacer de tal reacción termonu- 
clear. Un ingeniero del departamento de sanidad 
ambiental se encuentra observando las sutiles trazas 
de contaminación aparentemente presentes en una 
fotografía aérea de un río y se pregunta qué transtor- 
maciones podrá aplicar a sus datos para hacer avi- 
dente este problema de contaminación y poder tomar 

1po las decisiones correctas, Sobre la pantalla 
detelevisión acolor de aka resolución seven sorpren- 
dentes la primeras imágenes enviadas por el Voya- 
ger, en donde se aprecian por primera vez los 
acercamientos a mundos distantes de nuestro siste- 
ma solar. Las imágenes no obstante sondóbiles ylos 
científicos que las estudian dudan de algunos deta- 
llos, pues bien saben que el medio interplanetario 
pudo haber distorsionado signifcativamente lo que 











están viendo en la pantalla. Un grupo de geólogos 
está tratando de dilucidar el patrón de lineamientos 
observables en una imagen producida por el barre- 
dor multiespectral del Landsat-5, ya que la evolución 
tectónica de la zona está ligada a los fracturamientos 
y fallamientos que constituyen el patrón de lineamí 
tos, aunque si bien éstos se encuentran enmascara- 
dos en algunas partes del área cubierta pora imagen 
sateltaria. Un biomédico de un centro de salud 
comienza a comprender la evolución morfológica de 
Un grupede cólleinfctadas, elpatrón sapecialque 
de dichas células no se 
e e pp pol 
lo que el investigador biomédico requiere de más y 
mejores datos para eliminar las ambigúedades en su 
interpretación del fenómeno infeccioso. 


En todos los eventos arriba indicados, donde se 
describen diferentes situaciones en la investigación 
experimental, queda patente una problemática co- 
'mún: asaber, la caracterización matemática de datos 
denaturaleza bidimensional y el análisis interpreta: 
ción de los mismos. Para detallar esta problemática 
no olvidemos que en la generación y análisis de 
imágenes mutiespectrales, estamos hablando de 
hecho del diseño y la realización de un experimento 
bajo condiciones controladas y que en esto expari- 
mento estamos intentando medir un conjunto selecto 
de variables de estado que describen la situación en 
la que se encuentra un sísterna físico que deseamos 
estudiar, Además, el rango de medidas para dichas 
variables se encuentra también preestablecido, de tal 
forma, entonces, quelos datos que esperamos como 
resultado dela ejecución del experimento se encuon- 
tran bien definidos, en principio, en forma de una 
imagen muliespectral. 


Ahora bien, consideremos el siguiente esquema 
metodológico en el estucio de un sistema físico. Fig. 
51). 

Podemos apreciar en esto esquema, que en un 
proceso de retroalimentación, el comportamiento del 
sistema físico es estudiado a través del análisis de la 
imagen multespectral resultante del experimento, y 

















Sistema 
fisico(3D) sensor 

+ * + 
Comparación 


+ Escena 


t Modelo del 
sistema físico 


de datos bidimensionales 





Dispositivo 


Fenómeno de 
degradación 





Fig. 5.1 Esquema metodológico de estudio e interpretación de un sistema físico por medio del análisis 
multiespectrales. 
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que este análisis interpretado correctamente, sienta 
lasbases para generar un modelo del funcionamiento 
dedichosistemafísico. Comolacomparación de este 
modelo o simplificación del “mundo real” es factible 
por métodos cuantitativos, podemos entonces deter- 
minar el grado de adecuación que tiene el diseño de 
nuestro experimento para permitimos entender la 
realidad de nuestro mundo físico. No ahondaremos 
más aquí en este punto, pues esto constituye una 
pane fundamental del método científico aplicado alas 
ciencias experimentales y desde luego la Percepción 
Remota no escapa a este marco de referencia que 
constituye la metodología de la generación de cono- 
cimiento nuevo. Síppodemos decir, sin embargo, que 
el método cientifico establece claramente las bases 
para poder definir qué o cuáles escenas son necesa- 
rias capturar yregistrar y también nos permite estable- 
cerlas variables de estado y los rangos de medidaco- 
rrespondientes. La Percepción Remota es una meto- 
dología de observación a distancia de sistomas fs 
cos microscópicos o macroscópicos, por lo que es 
rosponsabllidad de cada investigador el definir su 
“mundo real" y las condiciones de investigación, así 
como las metas y objetivos correspondientes que so 
persiguen al establecer un experimento dado concer- 
niente a ese mundo real, Do aquí que nos ocupemos 
en la Porcepción Remota de la caracterización y 
análisis delcispositivo sensor, delas posibles fuentes 
de degradación y de la imagen multiespectral corres- 
pondierte. Aún cuando en el esquema del fgua 
estos aspectos se encuentran dilucidados, 
toda [laca an e imagen. 
Valde una escena dada, representa ho nada más la 
distribución espacial de las propiedades físicas pro- 
medio de ésta observadas a distancia, sino que 
también representa en forma implícita o explica la 
respuesta finita del dispositivo sensor y de los meca- 
nismos de interferencia, siempre presentes, a la 
medida de las variables co la escena. La imagen por 
tanto es una simplificación de la escena, limitada por 
el comportamiento delsensory porlos fenómenos de 
interferencia. 
mos, obtener 
interferencia o degradación, constituyen fenómenos 
fsicos bien defnicos que deben verse en realidad 
comosistemas físicos interaccionando, aveces estre- 
chamente, con el sístema que es de nuestro interés 
básico. Encientas ocasionesla degradación que sutre 
laimagen resultante no es de consideración; en otras 
'ocasiones.es detal magnitud que enmascara o distor- 
siona los datos producto del experimento, atal grado 
quela interpretación que se pretende hacer del siste- 
'mafísico setorna muy difícil .imposible. La; de 
de mecanismos de degradación en el proceso de 
captura deuna imagenno constituye unalimitación o 
Undefecto de las situaciones experimentales existen- 
tes en Percepción Remota, más bien es un hecho 
inherente a las ciencias experimentales en cuanto a 
los métodos e instrumentos de observación y en 
cuanto alcomportamiento de lanaturaleza misma. Es 
responsabilidad del investigador, sin embargo, el 
procurar evitaro minimizar cualquier tipo de degrada- 
ción en el arreglo experimental correspondiente. De 














loarriba expuesto debemos portanto dejar claro que 
alcaracterizar matemáticamenta una imagen multies- 
pectral hemos de plasmar también la respuesta fin 
del detector y las fuentes de degradación presentes 
Es por tanto la caracterización matemática de una 
imagen muliespectralel primer paso hacialainterpre- 
tación y modelado del comportamiento de nuestro 
sistema físico, 





Enla figura 5.1 se presenta un elemento complejo. 
dela Percepción Remota y es aquel que se refiere al 
análisis dela imagen muliespectral. En el capitulo 31 
introdujimos en forma conceptual los elementos más 
¡importantes para el análisis de una imagen digital yen 
se capítulo establecimos las bases para afirmar qu 
slanálsis de unaimagen se desglosa en dos grandes 
temas: el procesamiento de los datos multiespectra- 
les y el reconocimiento de patrones, de tal forma que 
podemos establecer la siguiente suma simbólica: 














bo 


Podemos, a su vez, decir que el procesamiento de 
Imágenes es el conjunto detransformaciones alespa- 
dio, donde se encuentran referidas originalmente los 
datos muliespectrales, con el objeto de poner en 
evidencia o realzar un cierto conjunto de patrones 
espaciales o espectrales. Es posible definir. cor 
mayorformalidad este conjunto detranstormaciones 
Para esto consideremos la figura 5.2, donde se 
resume el siguiente razonamiento: los datos que 
componen la imagen muiespectral original son gra- 
ficados en un espacio de tantas dimensiones como 
bandas tenga la imagen. 








Enla figura 5.2 se ha representado, por simplici 
dad, un espacio Euclidiano de tan sólo 3 dimensio- 
nes, pero en generales tendrán espacios «dimenai 
nales generados por sistemas de coordenadas ob! 
cuas, Los dominios %, + 0, están constituidos por 
jones conectadas simples y pueden verse tam- 
29m hipervolúmenes en espacios n:-dimensio- 
nales. Los operadores ', siendo acotados y finitos, 
realizan un mapeo de los datos del espacio de entr 
da (t.1.:) alespacio de salida(»,1.:, donde los 
datos están representados por las funciones, detal 
manera que en este nuevo espacio quedan en 
videncia o son reaizados cierta ciase de patrones. 
Desde luego, la forma de los operadores (%,) debe 
estar diseñada como para realzar precisamente 
aquellos pationes que son de nuestro interés 
Podemos decir aqui que el procesamiento digital de 
imágenes es aquel conjunto detécnicas matemáticas 
para el diseño y formulación de los operadores (0, 
necesarios en una O varias tareas específicas de 
realces de patrones, es decir que el procesamiento de 
imágenestiene como finalidad última elrealce de una 
ovarias clases de cobertura". Ahora podemos vercon 
mayor claridad, que en la serie de ejemplos dados al 
principio de este capítulo, los problemas comunes a 
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Fig. 5.2. Mapeo de los datos originales releridos al espacio ( : -.») y caracterizados por las funciones 
1, por medio de los operadores O, al espacio (xy.z) donde los datos están caracterizados por/as funciones 
9 Las funciones / y están definidas en los dominios D, y D, respectivamente. 








los quese enfrentanlos diferentes investigadores son 
lacaracterización de sus datos y.l realce de patrones. 

El procesamiento de imágenes, siendo un conjun- 
to de técnicas matemáticas y computacionales, pue- 
de ser dividido en tres secciones que son desglosa- 
das en razón del tipo de mapeo que ejecutan los 
operadoras (0), ostos tres grandes grupos se mues- 
tran acontinuación a laluz de lo que hemos discutido 
hasta ahora. 

La primera etapa pues en el análisis de datos expe- 
rimentals bi-imensionales consiste en la caracter- 
zación matemática de los mismos y una vez con esto 
'én lamano, se procede adar los primeros pasos para 
lanálisis que necesariamente está conectado con el 
tipo de datos con que se cuenta. Sólo los datos de 
carácter más general pueden ser realizados en cual- 
quiera de los tres tipos de realces mostrados en la 
figuras esto deberá quedar claro al finalizar 
la discusión del presente capítulo. Los realces unita- 
rios están constituidos por aquellos operadores que 
definen transformaciones unitarias o reversibles, es 
decir, los operadores correspondientes cumplen con 
un juego de condiciones de ortonormalidad que 
serán expuestas en el próximo capítulo. Por su natu- 








raleza, los realces unitarios no añaden ni quitan infor- 
mación al conjunto de los datos transformados, o en 
tras palabras, el mapeo que se esquematiza en la 
figura 52. puede ser hecho en ambas direccion: 
Los realces unitarios pueden verse también como 
tuna multirotación y elongación de los ejes coordena- 
dos del espacio original. En ciertas ocasiones los 
roalces unitarios se ven como transformaciones de 
tipo general para establecer criterios básicos | 
procesamiento ulterior, poro que puede decirse que, 
en este contexto, forman un grupo de preprocesa- 
miento de imágenes. Nosotros preferimos incluiralos 
realces unitarios como simplemente un grupo más en 
el procesamiento de 

A iferencia de los realces unitarios, los realces 
selectivos sí suprimen información, pues como su 
nombre lo indica el realce se realiza en forma select 
va,porlo que una clertaciase de patrones se pone en 
videncia a costa de otra clase, o bien a costa del 
resto de la imagen. En otras palabras, ciertos patro- 
nes son realzados y otros son suprimidos y como en 
este proceso una cierta cantidad de información es 
eliminada, resulta claro, entonces, que los realces 
selectivos no pueden serreversibles y desde luegono 
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Fig. 5.3 Primeros pasos para el análisis general de datos bi-dimensionales multiespectrales, 
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son unitarios. En esta clase de realces parte de los 
datos es separada de lanube original y bajo condicio- 
nes lineales las partes separadas pueden reunirse 
para recobrar la imagen primitiva. Sin embargo, en 
esta discusión suponemos que la parte separada es 
sliminada de hecho, por lo que la nube original es 
disminuida y por tanto la transformación resulta irre- 
versible, Losfítros de patrones espaciales y espectra- 
les son los exponentes más claros de los realces 
selectivos. Tal vez a primera vista podríamos decir 
que las reglas de segmentación expresadas básica- 
mente como una partición de los valores de los 
pixeles originales »,, tales que 


caen dentro del grupo de realces selectivos, ya que el 
valor del pixel de salida ”; está definido binariamente 
deacuerdo las funciones « delos pixeles de entrada 
ydeacuerdo aun umbral * de decisión. Sinembargo, 
aquellos pixeles de salida cuyo valor asignado sea 
uno no son en realidad suprimidos, sino que son ló- 
Sicamento separados delrestoyde acuerdo ala regia 
ú pueden ser eventualmente reunidos con el resto 
delospixelsoriginales. En este punto podemos hacer 
una afirmación muy importante: que todas las reglas 
detransformación puntuales, comola segmentación, 
son reversibles; pero que las regias de transforma: 
ción regionales son no-reversibles. Las transforma- 
ciones puntuales de hecho son combinaciones inea- 
les de los valores de los pixels en las diferentes 
bandas que componen la imagen de entrada y las 
transformaciones regionales se reñeren a operacio- 
nes sobre un cierto conjunto de pixels, normalmente 
conexos o vecinos, cuyo resultado produce un solo 
pixel de salida. Alas transformaciones regionales se 
lesllama también transformaciones de vecindad. Por 
supuesto que los fitrajes pueden expresarse como 
transformaciones de vecindad? las que generalmen- 
te se ejecutan banda por banda, aunque en ciertos 
casos puedan ser definidas sobre hiperregiones 
contenidas en el espacio generado por la imagen 
multiespectral original 


Enlosrealces de fidelidad se agrega informacióna 
laimagen original por medio de un modelo preesta- 
blecido y al cual se desea que se aproxime lo más 
posible dicha imagen. En este proceso sesuponene- 
cesariamente que en la captura y registro de la esce- 
na, la imagen correspondiente queda degradada por 
un mecanismofsico que se conocen principio y que 
se establece cuantitativamente como un modelo fsi- 
co-matemático. Se define entonces una métrica en el 
espacio de laimagen modelo y minimizando un error 
definido através de consideraciones estadísticas, se 
aproxima la imagen degradada lo más posible a la 
imagen modelo; esta aproximación puede realizarse 
¡en algunos casos enforma iterativassiempre y cuando 
la regla de minimización establezca un camino de 








convergencia que sea estable. En los realces de fide- 
lidadse procura, portanto, que la imagen degradada 
se parezca lo más felmente al modelo de imagen 
preestablecido y es precisamente este modelo el que 
inyecta información a la imagen degradada. Al acer- 
¡car nuestra imagen degradada al modelo ideal esta- 
mos otorgando una mayor calidad ala imagen trans- 
formada, pensamos implícitamente que la imagen 
resultante ha sido restaurada al estado en el cual 
estaría si no hubiese existido un mecanismo de de- 
'gradación, de aquí que a un realce de fidelidad se le 
llame también restauración de la imagen degradada. 
Es necesario dejar bien claro que la imagen no- 
degradada obviamente es desconocida pixel a pixel; 
lo que se conoce, entodocaso, sonsus propiedades 
estadísticas promedio ylo que es conocidocon cierta 
precisiónes el mecanismo de degradación plasmado 
en el mencionado modelo físico-matemático. Con 
esto queremos expresar que en la restauración de 
imágenes, lo que hacemos es definir un ideal de 
imagen, al cual deseamos que se asemeje lo más 
posible la imagen degradada. A este modelo o ideal 
deimagenle concedemos un valor mayor de acuerdo 
a las condiciones de nuestro problema. La imagen 
restaurada o no-degradada es obtenida entonces a 
partir de la imagen degradada, del modelo de degra- 
dación y delaregiade minimización de error. Quétipo 
de información es úti para nosotros y cuál no lo es, 
tene mucho que ver con el valor que asignemos a la 
imagen restaurada. El valor de la información está, 
pues, ligado a la calidad que pensamos tiene la 
imagen restaurada porque las nociones de valor y 
calidad son subjetivas al analizador de imágenes. Al 
hacer cualquier tipo de realce sobre una imagen 
multiespectral, estamos transformando los datos orí- 
ginales de acuerdo al esquema de la figura 5.2 
y los estamos enviando a un espacio donde en su 
"nueva representación les asignamos un mayor valor, 
Aunque suene paradójicamente extraño, esto es así 
aúnenlosrealces selectivos dondese pierdeintorma- 
ción, puesto que afin de cuentas estamos restringien- 
do voluntariamente nuestro universo de análisis a 
sólo una part de los datos originales mutiespectra- 
les y es precisamente a la nueva representación de 
esta porción de los datos a los que concedemos un 
mayor valor y alresto de los datos les damos poco o 
nulo valor porque no son de nuestro interés para un 
problemadeterminado. Entodos los casos del proce- 
samientode imágenes multiespectrales, se busca por 
tanto que los operadores (0, | sean tales que en el 
espacio de salida, la representación de los datos o 
parte de ellos tengan un valor mayor de acuerdo a un 
Problema dado y claro está no es posible diseñar un 
conjunto de operadores de propósito específico si no 
se han caracterizado adecuadamente los datos muk- 
espectrales. 


Ala luz de lo que se ha presentado hasta ahora, 
puede afirmarse que uno de los intereses básicos en 
el procesamiento de imágenes, se encuentra en la 
transformación conveniente de los datos multiespec- 
trales para relacionarlos a un marco de referencia 
adecuado auna interpretación deseada. Un segundo 





SELPER sec o —————————— o 13 


interés se encuentra en la extracción panicular de 
ciertosrasgos o características de una imagen, con el 
objetode apoyar un análisis humano oautomático de 
tales datos. En los dos casos mencionados, parte o 
todalaintormación inherentealaimagenesrealzada, 
inhibida, minímizada o simplemente eliminada con el 
objeto primordial de poner en evidencia uno o varios 
rasgos de la imagen determinados de antemano. El 
procesamiento de imágenes, así, se conviene de 
hecho en un ayudante inteligente para la intepreta- 
ciónya sea cualitativa o cuantitativa de la información 
contenida en la imagen. Hay que decir aquí que la 
interpretación cualitativa no necesariamente es subje- 
tiva, sino que, entodocaso, es un procesonoreferido 
auna eecala o patrón previamente definido. La inter- 
pretación cualitativa puede ser, por ejemplo, de tipo 
booleana, es deci, la decisión puede consistir en 
establecer la existencia o no de un patrón, o bien en 
lautiidadono de unaintormación dada. Silainterpre- 
tación es detipo cuantitativa, entonces puede referir. 
se ala existencia de uno o varios estados espaciales! 
de laimagen o ala morfología de una cierta clase de 
patrones En esteúlimotipo de interpretación se trata 
de decidir de acuerdo a una regla pre-establecida, 
sobre la presencia o ausencia de un patrón o una 
clase de patrones (clase de cobertura) y estimar o 
evaluar ciertos de los rasgos morfológicos de los 
patrones a través de un modelo matemático y una 
métrica. Los rasgos morfológicos de un patrón espa- 
cial o la manifestación espacial de un patrón espeo- 
tral, pueden referirse a la forma, textura, tamaño, 
perímetro, centroide, área, o parámetros geométricos 
de estructuras espaciales regulares. La decisión 
cuantitativa a la que nos referimos se conoce como 
reconocimiento de patrones y constituye el segundo 
gran aspecto del análisis de imágenes junto con el 
Procesamiento de las mismas. Én este trabajo no 
intentamos atacarel problema del reconocimiento de 
patrones por ser un tema extremadamente extenso y 
complejo. Enreferenciaal análisiscleimágenesahora 
¡podemos verlo de otra manera, por lo que diremos 
que este análisis consiste en la media morfológica 
de patrones, antecedida por la transformación geo- 
métrica de la imagen que los pone en evidencia. 
Entonces, el procesamiento de imágenes es aquel 
conjunto de transformaciones geométricas diseña- 
das para realzar uno o varios patrones, mientras que 
el reconocimiento de los mismos sólo puede hacerse 
através de los modelos matemáticos que caracteri- 
zan los diferentes aspectos morfológicos de los pa- 
trones, El reconocimiento de patrones sigue general- 
mente a la transformación geométrica, de tal forma 
quesia éstelo expresamos simbólicamente como ". 
y al procesamiento de imágenes como ',, entonces 
tenemos quel análisis de imágenes se puede expre- 
sar como: 





MáLisiA de Amiens = 0,10,(e 


es decir, el operador %. opera sobre la imagan 5 
después que el operador ',. y a esta sucesión de 


operaciones se refería la suma simbólica planteada 
con anterioridad en relación al análisis de imágenes 
(69. 
Elprimerproblema, pues, dentro de latemática del 
análisis de imágenes, se refiere a la comprensión de 
la estructura del sistema físico (fig. 5.1) que se desea 
estudiar. Esta comprensión implica finalmente la ela- 
boración de un modelo del sistema y que necesaria- 
mente tras consigo una simplificación o idealización 
de los parámetros o variables de estado que lo def- 
nen y, como ya se había mencionado anteriormente, 
el principio de la comprensión del sistema fisico 
estribaenlacaracterización matemática dela imagen 
multiespectal, conjuntamente con la respuesta finita 
del sensor remoto y los posibles mecanismos de 
degradación en la captura y registro de la escena de 
interés. Existen dos enfoques para la solución al 
problema de la caracterización o descripción mate- 
mática de imágenes muliespectrales: una descrip- 
cióndeterministay una descripción estadística detipo 
estocástico. Eltipo de descripción que se aplique ala 
imagen depende básicamente de su complejidad y 
del grado de simplificación empleado o deseado, En 
una descripción determinista se suponen conocidas 
o se conocen las relaciones de valor que generan 
¡cada pixel de la imagen y larelación con sus respec- 
tivos vecinos y, a partirde éstas se define una función 
matemática = = s(=.».1)donde se consideran enton- 
úcesrelaciones puntuales y; 'delaimagen. El 
parámetro. (iota) representa as diferentes bandas de 
la imagen y » está definida + (».1) < 9, De hecho 
podemos decir que al caracterizar la imagen, esta- 
mos caracterizando la escena correspondiente, con 
lo que para caracterizar al sistema físico debemos 
describir, entonces, todas las escenas relevantes a 
éste La función »11.".:) eslarelación que define 
la imagen correspondiente a una escena dada, es 
decir, la imagen es la representación matemática o 
'modelomatemático de una escena delsistema físico. 
Enlapráctica, sólo aquella sub-ciase deimágenes la- 
'madas gráficas pueden ser contoda certeza desciltas 
enforma determinista y en casos sencillos como imár 
¡genes binarias oimágenes tipo mapas, ladescripción 
Puede ser también la misma. En una descripción 
estadísticano se conocen, o no se emplean, las rela- 
'cones de valor que den una función definida y única 
entre los pixels de laimagen, sino que ésta se descr- 
be más bien a través de propiedades promedio, 
utilizando las funciones de distribución y densidad, 
funciones que desde luego son detipo probabilístico. 
Enuna descripción detipo estadístico no se conocen 
las relaciones de valor que generan cada pixel de la 
imagen; en otras palabras, a la hora de realizar las 
medidas sobre los CIV de la escena no se sabe con 
certeza cuál será el valor resultante de tales medidas, 
por lo que el valor que pudiera tener un pixel es 
manejado únicamente en términos probablísticos. 
Debe quedar claro quo lo quo os de orden probabilís- 
tico es el resultado de la medida sobre un CIV de la 
escena, puesto que ésta, por su complejidad, no se 
¡conoce y lo que sí es determinista esla respuesta del 
sensor, la cual es siempre la misma o al menos se 
¡conoce la manera en que responde bajo diferentes 
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condiciones. En el diseño y construcción de un sen- 
sor remoto, el estudio de la respuesta en diferentes 





que es 
¡sensor cuando un cierto tipo de predeterminada ra- 
diación incide sobre él. En una gráfica se conoce la 
función matemática que la caracteriza, por lo que se 
Puedesabercon probabilidad iguala uno cuál será.l 
valor de un pixel dado, es decir, se conoce larelación 
de valor que genera cada pixel de la gráfica. En » 
descripción estadíaticala variablo. == aí%.+ 

define como una variable aleatoria. En referencia aa 
figuraS.1 y alo expuesto hasta ahora, podemos decir 
'quelacaracterización delos diferentes aspectos enla 
aplura y registro de una imagen es como se muestra 
na siguiente tabla: 
























































Una imagen multiespectral se obtiene siempre a 
travós de un mecanismo de captura de información 
dela escena; los posibles mecanismos son muchos, 
pero pueden dividirse en directos (de contacto) o 
indirectos (Percepción Remota), dependiendo del 


método de captura. Sien este proceso de captura el 
valor de == síx.1,1) está bien determinado, el 
proceso es determinista; pero si estáregido por leyes. 
probabilíticas, entonces el proceso es estocástico. 
En cualquier caso dicho proceso puede verse como 
una operación realizada sobre la escena y de aquí 
que el sistema de captura sea entonces un operador 
que transforma la información del marco de referen- 
cía de la escena al de laimagen; o bien, visto de otra 
forma, la escena puede considerarse como una 
Imagen continuaqué transformada auna representa- 
ción digital queda registrada quasi-permanente. Esta 
consideración es muy imporanto ya quel transfer, 
imatemáticamenteig 


mación escena -imagen será. 
ocio 


ala transformación imagen « 
posible afirmar que una imagen multiespectral es 
siempre la representación espacial de la distribución 
de energía radiante proveniente de la escena para un 
conjunto de intervalos de energía, cada intervalo 
correspondiendo a una banda de la imagen mulies- 
pectral. Por energía radiante entendemos aquí, radia- 
ción corpuscular o radiación ondulatoria (acústica o 
electromagnética); es decir, nosreferimos a cualquier 
fenómeno detranspore de energía. En su forma más 
'General, la función que caracteriza dicha distribución 
aspacialde energía raianto depende delas coorde- 
nadas espaciales, del tiempo : y de un parámetro 
1. (Gota). En al caso de radiación ondulatoria, — re- 
presenta un intervalo de longitud de onda y enelcaso 
deradiación corpuscular, es un intervalo do energía 
, Sel haz de partículas incidentes al sistema físico, El 
parámetro, seportato una banda cualquiera ce una 
Imagen multiespectral. En sistemas modernos de 





extras o vinuales, no originalmente capturadas por el 
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Tabla 5.2. Constricciones y limitaciones 
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mentrno con Independencia obeervaciona, se tene 


al proceso de transformación de imágenes. 
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¡sensor remoto sino agregadas a la imagen captura- 
da. Así, por ejemplo, en estas bandas vinuales pue- 
den quedar almacenados datos geoquímicos, mag- 
néticos, gravimétricos o topográficos, en el caso de 
aplicaciones geofísicas, Para aplicaciones biomédi- 
cas, las bandas extras pueden contener datos quími- 
cos y para imágenes relacionadas con problemas 
industriales (radiografías y neutrografías), la distribu- 
ción de temperaturas superficiales bien puede que- 
daradicionada ala imagen original. Entonces escribi- 
remos la función que representa una imagen muties- 
pectral como. 


Bay E HOR SAY 4100 E 


Los valores que puede tomar la función g son 
finitos y dependen de la escala utlizada para repre- 
sentarlos. Los valores de las coordenadas (x.y) tam- 
bién están acotados y su intervalo de validez depende 
del tamaño espacial asociado a la imagen y del 
sistema de coordenadas empleado. En cambio el 
parámetro puede tomar en principio cualquier valor 
a parir de cero. Las consideraciones para eltiempo : 
son algo diferentes: se considera aquí que * es el 
momento en el cual ss captura la escena y para que 
esta consideración sea válida, se necesita que el 
tiempo que lleva capturar dicha escena sea mucho 
menor al tiempo que requiere el sistema fisico para 
axparimentar cualquier cambio. Con esta considera- 
ción estamos asumiendo enla práctica, que la captu- 
ra de una escena es “instantánea”, de otra forma 
tendríamos no-inealidades que no són convenientes 





sano también suponer que el sistema de captura de 
imágenes es no destructivo, sobre todo cuando se 
requiere de varias observaciones sobre la escena 
para obtener una imagen muliespectral, tantas ob- 
servaciones como bandas tenga la imagen. La fun- 
ción + puede ser vista como un campo escalar o 
vectorial. Sila imagen es muiespectral en relación al 


Fig.54.- Geometría de larelación del sensor 
con la escena. La cara intema s del sensor 
subliendealCIV cuya áreas S, sobreunalinea 
de barrido para un campo de vista (COW. 





parámetro +, es decir si + está definida para un 
conjunto de valores de :, elcampo es vectorial Pero 
si está definida para un solo valor de :, entonces el 
campo es escalar. La función + puede ser real o 
compleja dependiendo de si se mide o no lafase y la 
amplitud de una señal ondulatoria o de un haz de 
iones polarizados. 


52.- Proceso de Medida Sobre un CIV 


Enelproceso de muestreo para capturar y registrar 
una escena (fgs.11" . y5.4), uno de los primeros 
problemas que puede llegar a ser de consideración, 
es el provocado por la apertura finita (»)del detector 
de un sistema de captura de datos bi-dimensionales 
que, a través de un mecanismo de barrido (punto a 
punto, ineal o matricial) produce una imagen digital. 

El ángulo sólido finito definido por la apertura finitade 
este detector, introduce una frecuencia de corte que 
debe ser tomada en cuenta cuando se establece la 
frecuencia de muestreo (Teoremade Muestreo). Esto 
es importante, ya que dicha frecuencia de corte pue- 
deservistacomo elresultado de un ftraje pasa bajas, 
el que en algunos casos puede ser utlizado con 
ventaja para eliminar patrones oruido de altafrecuen- 
cía que no son de interés para el observador. Haga- 
'moseste análisis en dos pasos; en primera instancia 
consideremos la frecuencia de corte introducida por 
la apertura finita del detector, cuando ésta no es muy 
pequeña comparada con el área S del CIV, para 
posteriormente obtener la frecuencia de corte produ- 
ida por el valor finto del CIV. Además de esto, el 
¡campo instantáneo de vista, el que está en relación 
directa con la resolución espacial del sistema de 
captura, depende de la geometría particular escena/ 
detector. Es por esto que en algunas ocasiones, al 
CiVselellama elemento de resolución. Sea entonces 
selárea interna del detector sobre el cual se proyecta 
elcampo instantáneo de vistay sea (u,.)elsistema de 
coordenadas cartesianas a las que s es referida. A 
medida que evoluciona el proceso de muestreo a lo 
largo de una línea de barrido, el detector capta sobre 
su cara intema, una distribución de intensidad de 
radiación lu. .:).De esta manera y como una res- 
puesta del mecanismo de detección, se establece la 
salida del detector una distribución de valores (nor- 
malmente votajes) de intensidad dependiente del 
tiempo, proporcional a la cantidad de radiación que 
incide por unidad de tiempo, por unidad de ángulo 
sólido y que está dada por 
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y que generalmente se mide en mw/s/sterad. A 
medida que pasa eltiempo se produce una secuencia 
de valores, - X*,Juno para cada CV de la lín 
barrido. La salida del detector normalmente es dis- 
cretizada por medio de un convertidor análogo/digi- 
tal, a una escala de valores aqui discretos con un 
rango de! » 2 - 1, donde » se escoge entre 5 y 8. 
Estos valores discretos son los que constituyen los 
pixels de la nea de barrido en cuestión. 
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Seaahoraylaabcisa al origen de lalínea de barrido 
yxlaordenadaal origen del CIV. Sea también (g.5.4) 
U paralela a + y 1-17 el intervalo cerrado que 
cubres sobrela dirección u. Entonces para un barrido 
lineal se tiene x - k:,conde k es una constante. 
Algunos sistemas de captura (Landsat) tienen barr: 
idos no-ineales pues utlizan un espejo osclatorio; 
paraóstos x= kses(w), donde wos la frecuencia de 
oscilación. No restamos generalidad al presente 
desarrollo cuando consideramos un barido lineal, 
puesto que sl mismo racionamiento puedo ser hecho 
paraotras relaciones entre» » *siempre ycuando se 
tonga la función correspondiente. Por tanto la distri- 
bucióndle intensidades h(+)puede servistacomo una 
función de x 














donde !' es la distribución de intensidades para un 
barrido lineal, Es decir, que se ha eliminado la dopon- 
dencia on eltiempo para proceder a posicionarse en 
una xcen la línea de barrido y sumar así sobre u para 
sl ároa que cubre la cara interna del detector, Con el 
cambio de variable w" = y + x/% y escriblondo expli 
citamente los límites de Integración para u!, se tiene 








peta 
w0 +) JE 





ya que du» du puesto que x se fa cada vez que so 
Suma sobre u. Detal forma que si se denomina, ala 
suma de intensidados sobre la coordenada y, 
entonces la distribución Ax) se escribe como 
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Hemos hecho esta simplificación pues suponemos 
que el barrido sobre x es independiente al barrido 
sobrey, además de queelanchodela líneadebarrido 
se considera ser mucho menor al tamaño de la 
imagen. Laoxpresión (5.8) nosindicaqueeldetector 
'esuna ventana que se mueve alolargo del eje x y que 
(4) esla distribución de promedios de intensidad de 
radiación sobre el mismo eje. El efecto que produce 
esta integración sobre x es el de un fitraje de paso 
bajo, con una frecuencia de corte que está explícita- 
mente relacionada con elrango. [-1/2,1/21 “elque 
se encuentra definido por la apertura del detector. 











Sea ahora ¡() latranstormada de Fourier de 








donde 1 es un índice mudo de suma, Consideremos 
también la transformada de Fourier de 41) 





Aa) 1/0 rx 
A 
60 
y con el cambio de variable —1' =1-x/k válido 


desde luego para el barrido lineal, se tiene 


» 1 
Ko) TRA) alex 





y tooscriblendo ésta úhima exprosión, nos queda la 
siguiente relación 
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¡Como ahora transformada de Fourier 
de! entonces por el teorema de Shifing tenemos lo 
siguiente 


02, 


A 
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(513) 


Lo que quiera decir que después de la integración, 
la expresión para la transformada de Fourior dl la 
distribución de intensidades queda finalmente como 








(a1/2K)/u = Tía) 1 Hne(ot/2) 





Tanto como  seencuentranen el espacio 
de frecuencias lo que permite visualizarla frecuencia 
de corte del fitro pasa bajas introducido por la 
apertura finita del detector. Un análisis de la función 
pinciol/2k) — muestraquepara, lol = 2/1, 
la distribu-ción A(s) es atenuada y prácticamente 
reducida a un valortal que puede considerarse como 
ero, puesto que el 22 lóbulo de la función sinc(x1/2K) 

del orden de un vigésimo con respecto al 
1*lóbulo, Estoindica entonces que para una apertura 
finita y no despreciable con respecto al CIV, se 
introduce una frecuencia de cone que es 
inversamente proporcional al tamaño de la apertura 
deldetector. De aquí que tengamos un primer princk- 
plo de conservación que nos dice que la información 
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en 1(x) no puede exceder a la de la escena. Otro 
punto importante en relación a estas considera- 
ciones, en el proceso de transformar una escena a 
na imagen digital, es el de saber con que frecuencia 
durante el barrido, la función A(x) va a ser evaluada; 
esto será tratado en el Teorema de Muestreo. ES 
esencial conocer esta frecuencia de muestreo, pues 
ella determina directamente el intervalo de tiempo 
óptimo porel cual se divide eltiempo total de barrido. 
Desde luego la frecuencia de muestreo no puede 
exceder a la frecuencia de corte -. derivada 
anteriormente y, además, esta frecuencia de mues- 
treo nopuede sertal que impida al detector permane- 
cereltiempo mínimo de residencia necesarioen cada 
CIV de la escena. De otra manera la relación señal/ 
ruido no seríala adecuada para el valor digitizado de 
AG). En otras palabras podemos decir que la 
frecuencia de cone, la de muestreo yla respuesta del 
detector son parámetros de suma importancia en 
todo diseño de un sistema de captura de imágenes. 


El segundo punto importante que vamos a consi- 
deraran relación ala geometría: sensor/escana, es el 
queserefere ala frecuenciadecorte ..'* introducida 
PY llama ato del CV; donde en general -- 
ya que s > » (19.54). El razonamiento es muy 
similar al arriba desarrollado, ya que el sensor está 
midiendo las propiedades promedio de todos los 
objetos que se encuentran en el CIV enfocado por el 
sensor. En otras palabras, el sensor integra toda la 
radiación, por unidad de tiempo, por unidad de 
ángulo sólido, que proviene del CIV en cuestión, 
siendo esta integración de radiación a la que nos 
referimos como las propiedades promedio. De esta 
forma se introduce una frecuencia de corte cada vez 
mayor a medida que se reduce el tamaño del CIV, 
puesto que los ocuparán objetos de menor tamaño. 
Con este razonamiento podemos afirmar sin 
necesidad de calcular la integración de raciación 
sobreelCIV. quela frecuencia de corte está dada por: 








5.15) 


donde L es el lado de un CIV cuadrado. Para CIV* 
rectangulares se tendrán recuencias de conte diferen- 
tes para las direcciones x,y. Es claro que en el diseño 
deun sistema de captura de imágenes lo que importa 
es saber cual es la frecuencia de corte menor: la del 
sensor o la del CIV. En la investigación experimental 
'es muy posible que lo que defina finalmente eltamaño 
de la apertura del detector y del CIV, sea el conjunto 
derestricciones de índole práctico como: la distancia 
del sensor a la escena, la respuesta del sensor, el 
movimiento relativo del sensor y la escena y, el volu- 
men de datos razonablemente manejables por com- 
putadora. Sólo en contados casos setiene suficiente 
libertad para decidirla frecuencia de corte que mejor 
convenga a los objetivos de un experimento dado. 
Hemos analizado hasta aquílos aspectos más impor- 


tantes de la mecición de las propiedades promedio 
de cada CIV de la escena. 


Supongamos ahora que tenemos un sistema de 
barrido para discretizar toda la escena, veamos en- 
tonces cuantos CIV debemos medir, en que tipo de 
arreglo espacial deben estar estos elementos sobrela 
escena y que condiciones se deben cumplirpara que 
slconjunto de medidas represente adecuadamentea 
la escena. 


53.- Teorema de Muestreo 


Un segundo aspecto de gran importancia en la 
captura de una escena, es como ya se había men- 
cionado antes, el de la frecuencia de muestreo para 
discretizar el promedio de intensidades A(). Esta 
frecuencia de muestreo debe ser tal que permita la 
reconstrucción del detalle espacial más fino de la 
escena. Por otrolado es esencial muestreara función 
Acx)a la frecuencia Óptima, pues un sobremuestreo 
produciría un exceso de datos innecesarios que ha- 
rían más costoso el trabajo computacional. Además 
siel muestreo es subóptimo se perderán detalles de 
la escena que pueden ser importantes en la nterpre- 
tación de un fenómeno dado. Es importante mencio- 
nar que ciertas aplicaciones del análisis de imágenes, 
como la prospección geofísica a partir de imágenes 
sateltarias, no requieren de la preservación de los 
detalles más finos de la escena, más aún éstos son 
innecesarios dada la escala regional requerida en 
dicha prospección. En otras aplicaciones como la 
Biomedicina se necestará de todo detalle posible, 
La frecuencia de muestreo mencionada la determi- 
na el Teorema de Muestreo” que se enuncia como 
sigue: “Una función 1.1) cuya transformada de 
Fourier es igual a cero para, > 
está univocamente determinaca por los valores 
tomados en puntos uniformemente espaciados en el 
plano (x.y). siempre y cuando los espacios entre los 
Puntos satisfagan las squietes condiciones: as 
1/2. En términos físicos, el período 
de muéstreo debe serigual o menor que la mitad del 
tamaño del más fino detalle dentro de la escena. Para 
clarificar esto un poco más recordemos el concepto 
de frecuencia espacial desarrollado en el capítulo 3!, 
en donde se explicó la idea de textura de una imagen 
digital. La digtización espacialde una escena implica 
necesariamente su cuantización espectral, la que 
consiste en establecer un intervalo del espectro 
electromagnético o de energía, en el que se hacentas 
medidas correspondientes y en definir una escala de 
Niveles discretos contra la que se comparan las 
intensidades A(+)de radiación provenientes de cada 
elemento (CIV) de la escena. Al proceso de 
discretización espacial, es decir a la división de la 
escena en CIV regulares y no separados le llamare- 
mos digitización. Al proceso de discretización espec- 
tral le Namaremos cuantización (también se le lama 
discretización radiométrica), Ambos son formalmen- 
te independientes, pero tienen mucho que ver en la 
representación adecuada de la escena. Entonces el 
proceso de discretización espacial tiene como resul- 
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tado una matriz 





TA 


Lea 


a 








para una banda + y para (¡) enteros. 
Nótese que el tiempo no aparece 
puesto que consideramos que lacaptura dela escena 
es, para todo orden práctico, instantánea. Los ele- 
mentos de la matriz (5.16) tienen las 
características: 1). El elemento f(j), de la imagen 
digital representa la subárea $x0y — delaescena 
0 imagen continua; 2).- Existe una relación univoca 
entre cada elemento de la imagen digital y cada CIV 
dela escena; 9). Elvalo del elemento 1) relerdo 
escala do niveles discretos, es proporcional alas 
características espectrales promedio de las medidas 
hechas sobre el CIV de la escena; 
4) El número de elementos de la matriz (5.16) para 
na escena dada, es decos valores do N y M, Jato 
conla escala de niveles discretos, definelaresolución 
de la imagen digital, os decir el grado de detalla 
disconible dela escena; 5). El número de elementos 
(NM), la escala de cuantización, la relación señal/ 
uldo, los factores de contraste! y la agudeza de los 
bordes de lasclases espectrales, definen en conjunto 
y en una relación compleja lo que se conoce como 
calidad de la imagen digital. Como en este concepto 
de calidad intervienen los mecanismos de la 
Psicovisión y laapreciación humana y portanto dicho 
Concepto cas en el campo de la subjetividad. Así, por 
ejemplo, una pintura de esto subrealista, con pocos 
elementos (N,M) y poca agudeza en los bordes, le 
asignaremos, a pesar de todo un alto de grado de 
calidad. O bien una fotografía borrosa de un auto en 
movimiento, sería probablemente elegida para una 
exposición de arte a pesar de no satisfacer los 
criterios de calidad arriba mencionados. En otras 
palabras, si bien la calidad de una imagen está en 
principio, definida por los parámetros arriba 
indicados, el modelo que utilizamos par 
'combinaros y producirun índice de calitad, depende 
de la aplicación panicular que deseamos darle a 
nuestra imagen digital. 











Elproceso de muestreo odigitización puede visua- 
lizarse con la ayuda de la sucesión de figuras 5.58 - 
55k. Por simplicidad y para entender primero en 
lorma cualitativa dicho muestreo, consideraremos 
que hacemos un corte, a lo largo de una dirección 
paralela al eje x, de la superficie que representa los 
valores de la escena, obteniendo de esta forma la 
función £(x) mostrada en la figura 55.2. En otras 





palabras es como sí esta !(x)fuese la sucesión de 
valores a lo largo de una línea de barrido, Como lo 
supone el Teorema de Muestreo, la escena está 
limitada en frecuencia, lo que quiere decir que la 
respectivatransiormada de Fourierdefla): — Fw), 
será poco más o menos como se aprecia en la igu 
5.5. Para obtener una versión muestreada de t(x 
basta con multiplicar los valores de una función de 
muestreo (xi (6ig.5.5.c), por los valores 
correspondientes de:»). La función sí consiste de 
untren deimpulsos, infinitamente estrechos, de altura 
constante ycon un espaciamiento ¿x entre ellos. Esto 
desde luego en la práctica no es posible puesto que 
el CiVno puede hacerse infinitamente pequeño, nía 
tespuesta del sensor puede ser infinitamonto 
eficiente. La transformada de Fourier de s(x) (fig. 
5.5.0), estambién untren de impulsos, pero con una 
distancia entre ellos de 1/3 
frecuencias. Entonc 
Multiplicación de »1x) con £(x) (Kg. 5.5c), en ur 
función discreta que representa ala función continua 
£10.La posiblidad de recuperar integramente esta 
función continua, a panir de la función discreta 
aíx)Híx), depende de la capacidad para poder 
Interpolar correctamente entre los valores discretos 
des(x!(s y recuperar asíla función continuat(»).Esta 
capacidad está dada univocamente por la distancia 
is, 0 equivalentemonto por la fracuencia de 
muestreo, La transformada de Fourier de la versión 
loca asco comcución Tasa). (ho 
554), la que por la naturaleza de la convolución” as 
sn realidad una sucesión de copias de, ilu) — que 
Pueden estar o no traslapadas dependiendo de la 
relación que guarde la frecuencia de muestreo 1/0: 
con la frecuencia de corte, Esclaroquesi 
la frecuencia de muestreo es relativamente alta. (fig. 
5.59), las copias de. 1) aparecerán separadas 
sintrasiape. Es claro también que la condición de no- 
tanapo qa tapan que, 0 5 l/s oben 
a Es decir que para poder 
eoorrlafunción conato pedi dos lores 
'muestreados, la frecuencia de muestreo debe ser 
'mayoro igualque el doble de la máxima frecuencia (o 
detalle más fino) contenida en la escena alo largo dal 
eje x. Un razonamiento similar al anteriormento ex- 
Puesto, puede hacerse para los valores de la escena 
alolargo del eje y, siempre y cuando el muestreo alo 
largo de un eje sea independiente del muestreo a lo 
largodelotro eje. Sila anterior condición de muestreo 
se cumple, entonces con una función escalón de 
valoradecuado (1.5.5), se pueden eliminar. todas 
lascopias de *í=) menos una (fg.5:5.k). Lo que 
algebraicamente se escribe como. 
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de donde obviamente se puede obtener la función 
¡continua £(x) aplicando la transformada inversa de 
Fouriera 5) El razonamiento cualitativo 
de las figuras 5.5a-k puede ser fácilmente extendido 
a dos dimensiones ya que en lugar de gráficas de 
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Fig. 5.5.- Diagrama cualitativo de los varios aspectos que componen al Teorema de Muestreo. 
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curvas se pueden utlizar gráficas de superficies. 


Varias limitaciones de orden práctico existen sin 
embargo al Teorema de Muestreo, que deben ser 
tomadas en consideración en el diseño de un sistema 
'decaptura deinformación bidimensional, Las restric- 
ciones en la vida real pueden resumirse como sigue: 
1). La función de muestreo «!») nunca es un tren de 
impulsos infinitamente estrechos; generalmente es 
una serie de gausslanas cuya anchura depende del 
sor remoto; 2). La función 1 x no está limitada en 
3).-La función 
160 está definida para ==: 0 < Todo 
sto hace que en realidad la frecuencia de muestreo 
tenga que ser un poco mayor a la indicada por ol 
Teorema de Muestreo, tomando en cuenta desde 
luego que la frecuencia de muestreo no puede ser 
mayoralafrecuenciadecorte impuesta poreltamaño 
finito del CIV y por la apertura del detector. 

Porotro lado la validez de la representación discro- 
tadelaescena (xy) depende de sis posible recons- 
truir la función original a partir de sus mu 
el objeto de obtener una demostración 
del Tooromade Muestreo, introducimos ahor 
de laticas de muestreo ya que hasta ahora sólo 
hemos hablado del tema en forma intuitiva, al supo- 

¡plcitamento una rela regular con base en la 
cual se realizan las medidas acerca de la escena. 
También hemos supuesto por comodidad que los 
CIV son rec y cuadrados, contiguos y sin 

¡pe No siempro es este el caso puesto que 
existen ocasiones donde el muestreo es nohomogé- 
neo, los CIV no sonrectangulares y existo un traslapo 
entre ellos, En cualquier caso el arreglo que guardan 
los CIV es siempre isotrópico pues de otra manera 



















































iboratorio, los. 
CIV son cuadrados o rectangulares y justamento sin 
traslapo, Ensistemas de captura a bordo de vehículos 
orbitales se dan variostipos de latices. Por ejemplo en 
los satóltes Landsat y Spot los CIV van incrementan» 
dosu tamaño a medida que se alejan del eje longitu- 
¿inal de la escena y esto se debe a la curvatura de la 
superficie terrestre, El barredor muliespectral del 
Skylab tanía un barrido detipocónico dexalformaque 
losCIV estaban formados porsectores circulares que 
incrementaban su tamaño a medida que se alejaban 
del eje de barrido, En ningún caso se tienen redes 
irregulares de CIV. 














Las ideas anteriores van a quedar formalmente 
establecidas por la introducción de las latices de 
muestreo, lo que además constituye una herramienta 
matemática muy útil para entender el Teorema de 
Muestreo bidimensional. Estas laticos se definen* en 
sl espacio de coordenadas y en el de frecuencias y 
¡son mutuamente recíprocas desde el punto de vista 
vectorial. Una latice es una estructura lógica que nos 
define la posición y distribución de los CIV sobre la 
escena. Lalatice es una rejila virtual que se sobreim- 
pone la escena y con la cual queda determinado el 
tamaño de un CIV de acuerdo a la posición relativa 
que guarda con respecto a los demás. Esta ltice se 
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Fig, 56. Latice regular de muestreo (a) generada 
poros vectores r.yr,y larespectiva late recípro 
Sa (b) definida por w, y Wo. 











¡genera por medio de una base vectorial bi-dimensio- 
nal, entonces sñan", ) dos vectores en el espacio 
de coordenadas (9.5.6) ortogonal.) Los sisto- 
mas de coordenadas mostrados en la figura 5.6 no 
necesariamente tienen que ser ortogonales pero es 
lo sean; lo que es posible es que 
jean no-euclidlanos, es decir poctemos tener sisto- 
mas de coordenadas oblicuos, A los vectores! .+2) 
les llamaremos base vectorial y constituyen la estruc» 
tura matemática con la cual generamos la latice y de 
aquíque óstala podamos ercomo un arreglo regular 
de puntos en el plano! 1). Estos puntos están del 
idos por los vectores de posición 1, dados por 











A 





(5,18) 


com OA 2. *1i para cuando 
slorigon de coordenadas se coloca en ol contro de la 
escena. En general los coeficientes (m,n) dependen 
der, Y Ez , pero aquí, por simplicidad en el 
razonamiento se suponen constantes, En el presento 
trabajo suponemos que el origen de coordenadas se 
"encuentra en la parte inferior izquierda por lo que m y 
1 los tomamos como números naturales, La base 
vectorial», v,) delalatice recíprocasn el espacio de 
frecuencias está univocamente determinada por los 
vectores base originales, cuando se supone un 
isomorfismo entre los dos espacios. La relación es la 
siguiente 

















El producto escalar dela relación (5.19) nos indica 
quelos vectores originales son mutuamente perpen- 
diculares alos vectores recíprocos, 
L «Por otro lado, de la relación (5. 
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donde es el ángulo que forman los vectores r, yr. 
La base vectorial recíproca, análogamente a la baso 
original, genera un arreglo periódico de puntos en el 
determinadas por el vector w,,, dado por 


Mea 





La latice del espacio de coordenadas constituye 
finalmente la manera como el detector va a ir midien- 
do las propiedades espectrales de la escena. O bien, 
la latice define la distribución de lugares (CIV) donde 
el sensor remoto medirá, por unidad de tiempo, la 
cantidad de energía proveniente de tal escena. 

Sea ahora una función bidimensional continua 
xy) que representa una escena cualquiera y sean 
£(i,j) las muestras correspondientes tomadas en las 
posiciones (ij) de la escena, donde el número de 
muestras define el tamaño (M,N) de la imagen digital 
correspondiente. Se supone además que las mues- 
tras fueron tomadas para un valor determinado del 
parámetro : y que el tiempo de captura es muy 
pequeño comparado con cualquier cambio que la 
escena pudiera tener. Por comodidad, esas funcio- 
nes serán escritas de ahora en adelante como 11) y 
1(r) respectivamente, Donde res el vector de posición 

¡punto (x.y) y 1, indica que es el vector de posición 
para valores distretos del plano x-y. Una manera 
conveniente de saber si 1(r) puede ser reconstruida a 
partir de las muestras (1,) es interpolando en forma 
Correcta entre ellas, sur 'claro está que dispo- 
'nemos de un número arbitrario de medidas sobre la 
escena, Más aún a partir de aqui se puede determinar 
matemáticamente la condición de muestreo que fue 
visualizada (fig. 5.5) anteriormente en forma gráfica. 
Si1(r) puede ser reconstruida a partir de f(r) entonces 
estamos implicando que la imagen es una 
representación válida de la escena, derivando al 
mismo tiempo las condiciones de validez. Sea anton- 
¡ces g(r) una función de interpolación, la que convolu- 
cionada con ((r,) genera la función original que repre- 
senta la escena. 








OS 





Es decir que intentamos recobrar la función 1(), 
multiplicando el valor delas muestras poruna función 
de peso adecuada y haciendo notar que todas las 
muestras contribuyen cada vez para generar un sólo 
punto de la función 1(). Además de suponer que 
existe g(1), asumimos también que existen las trans» 
formadas de Fourier de fy 9. Otra manera de vera la 
relación (5.22) es determinando las condiciones bajo 
las cuales la función f() puede ser expandida, en 
series, por medio de una relación lineal de sus mues- 
tras y con una función de peso gl). Ahora bien, por la 
propiedad de Shiing dela función * de Dírac,se 


reparte pue 


cm 


Substituyendo ahora esta ecuación en la expresión 
(522), e inviniendo el orden de las sumatorias y la 
operación de integración, dado que fy g son funcio- 
nes bien comportadas. La función original se expresa 





La función 
dirección x y con periodo 1, en la dirección 
quiere decir que puede ser expandida en 
Fourier 





> 7 E ajayesolanite! +9 dT,T,) 
Pra AT, 


donde los vectores *' 7 %,, pertenecient 
espacios de coordenadas y de frecuencias respecti- 
vamente, están dados por 








5.26) 
Lasvariables (xy) son continuas y las (4,1) discro- 


tas. Ahora bien, los coeficientes 1, y », de la expan- 
sión (5.25) están dados por 
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haciendo notar que %1 "0% Y Tay Td, 

hemos supuesto que el muestreo es independiente 
de una u otra dirección, coma ya lo habíamos 
anotado anteriormente. Aquellos términos con m,n? 
1 quedan fuera del rango de integración, por lo que 
sólo los términos con mn = O de la suma 
contribuyen a ésta, de tal forma que debido a la 
propiedad de Shiting llegamos al siguiente resultado 


dy VI, Ema 
62 
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Substtuimosahoraelresultado de(5.28) ylaexpre- 
sión (5.25) en (5.24) para tener una relación ciara de 
la expansión lineal de f() en términos de gí1), así 
entonces 





Al analizar de cerca esta Úlima relación podemos 
apreciar que se trata de la doble sumatoria de la 
convolución de las siguientes funciones 





de tal forma que por el teorema de la convolución, 
establecido como una de las propiedades de latrans- 
formada de Fourier, que dice que la transformada de 
Fourier de f(1) es igual al producto delas transforma- 
das de Fourier de las funciones involucradas en la 
convolución. Sin embargo, hay que tomar en cuenta 
que la primera de las funciones involucradas en 
(5.30), contiene un cambio de fase dado por la expo- 
únencial, de manera que en el espacio de frecuencias 
esto equivale aun corrimientoenta frecuencia. Enton- 
ces de acuerdo a las exprosionos (5-29) y (5.30), 
tenemos que la transformada de Fourier de la función 
que representa a la escena es 


a 


Es necesario estudiar ahora con detalle esta útima. 
relación pues de ella obtendremos las considera- 
ciones para la frecuencia Óptima de muestreo, así 
como las limitaciones respectivas. En primer lugar 
únolemos que T,, spresenta un área en el espacio de 
coordenadas cón un equivalente en el espacio de 
frecuencias, ya que se trata de la frecuencia de mues- 
treo en ambas direcciones: x.y. También hay que 
dejar claro que los puntos dados por », forman la 
latico recíproca de los puntos de muestreo que a su 
vez forman la ltice original. Refrámonos ahora a la 
figura 5.7, en dondese aprecian varias copias de la 
transformada de Fourier de “(r) para un conjunto 
selecto de valores de «.1, Las copias no se encuen- 
tran traslapadas porque la frecuencia de muestreo 
excede al valor óptimo explicado en la figura 5.5 y, 
dado que la:(r'se supone limitada en frecuencia, las 
mencionadas copias aparecen con simetría alrede- 
orde los ejes =,: »,. En forma análoga al caso uni- 
dimensional debemos encontrar la frecuencia de 
muestreo para que las '(»-x, no setraslapen para 
diferentes valores de k y 1, esto implica que es 
necesario escoger adecuadaente a w, y aw. los que 
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finalmente determinan r, y r,, los que finalmente 
determinan r, y1,. Para lá situación de exactamente 
noraslape, podemos escoger a 6(w) de tal forma 
que cubra una área igual a la que ocupa Fw) y cuyo 
valor sea el de, Para aquellas regiones no cubier- 
tas por tw) osuscopias, ((«) puedetomar cualquier 
valor, porlo que esta función de interpolación no esta 
univocamente determinada. Sin embargo, esto lo 
podemos tomar con ventaja y seleccionar aw) 
¡como un paralelepípedo para (NG. 5.7) que encierre 
justamente a *(w) y cuya base sea precisamente !,,. 
Es muy importante anotar que esto es válido para 
campos aleatorios homogéneos y limitados en tre- 
cuencia, ya que de otra manera hay que buscar una 









to q 








Fig. 57. Copias no traslapadas de la Trans- 
ormada de Fourier F(w=w,) de la función 1(), la 
que se supone limitada en frecuencia. Nótese 
que tanto (4) como su Transformada de Fourier 
son superficies y aquí se muestra tan sólo una 
vista superior de éstas. 









estrategia de minimización de error? para la recons- 
trucción det). Entoncesen formasimilar al caso uni. 
«imensional se tene que 1», a la que so la conoce 


¡como función fitro está dada por 


1, 
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DON Iv y Recordemos que el 
"espaciamiento correspondiente de los puntos mues- 
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treados en el espacio de coordenadas se obtiene a 
partir de la relación entre el período de muestras y la 
máxima frecuencia contenida en la escena, a la que 
se le conoce también como frecuencia de corte, por 
tanto. 





1, MUS Ey dy 





15m) 


donde bx == y br= 0 ¿Obien E Bxrie dy, por 
lo que con esto queda bien definida la atico de 
muestreo, Regresemos ahora ala relación (5.32) y 
escojamos de entre la infinidad de posibilidades, la 
más simple para 5(w), es decir aquella donde ésta 
función filtro está dada por un paralelepípedo (fig. 5.7) 
contados /-: y ¿»; y altura G. Por lo que sitomamos 
la transformada "de Fouier inversa de “lw 
obtenemos la siguiente expresión 





De esta última relación es claro que si escogemos 
al periodo de muestreo tal que se satisfaga 
precisamente la condición de no-traslape, entonces 
si 1-2: ysi 1,=1/% la función de 
interpolación se escribe como 





que es precisamente la transformada de Fourier del 
mencionado paralelepípedo. El razonamiento puede 
serhecho ale inversa, es decir, preguntaros cuando 
debe valer *., para ques!) sea una función sine bi- 
dimensional y de ahi obtener la condición de no- 
traslape. En el mismo orden de ideas podemos ver 
ácimento ques submuestreamos, es decis 

y 1, > 1/2 entonces no habrá manera de 
"encontar ná Gu) que sausega la relación (532) 
Por otro lado si sacamos más muestras de las 
necesarias la olección de G(w) puede seguirsiendo!a 
del paralelepipedo, Es necesario puntualizar que el 
desarrollo anterior para el teorema de muestreo bi- 
dimensional se refiere a una imagen de una sola 
banda o bien a imágenes multiespectrales descritas 
por campos aleatorios homogéneos (ver descripción 
estocástica y limitados en frecuencia, Esto se debe a 
que las condiciones de muestreo no pueden 
depender de las propiedades estadísticas de cada 








bandadelaimagen muliespectral yesto sólo selogra 
suponiendo. válidas las condiciones. arriba 
mencionadas. 


Se han introducido hasta aquí una serie de aspeo- 
tos cuantitativos relacionados a la caracterización 
matemática de una imagen y al análisis general de la 
misma. Se introdujeron condiciones locales de medti- 
da, es decir, estudiamos el proceso de captura de 
información sobre cada CIV y posteriormente vimos 
los aspectos regionales de captura de una escena al 
introducir el teorema de muestreo. Vamos ahora a 
resumir todos estos conceptos en un modelo mate- 
máticoque describelatransformación de una escena 
aunaimagen digital. A esta transformación la pode- 
mos entender también como la definición de varios 
operadores de captura de información quese aplican 
bajociertascondiciones y limitaciones.en el momento 
de estudiar una escena determinada. Es importante 
puntualizar que, dado que una escena puede ser vista 
como una imagen continua, el modelo matemático 
que obtengamos será de validez general. Debido a 
esto, escribiremos este modelo en varias versiones 
aplicables a diferentes situaciones donde estén invo- 
Iucradas imágenes continuas y discretas. 








5,4.- Operadores en Captura do una Escena 


La captura de una escena en una imagen tiene 
básicamente surazón de ser enlanecesidad dereferir 
cierta información de una escena a un nuevo marco 
de referencia a partir del cual, el tratamiento, extrac- 
ción ymodelado de patronesdde ésta, sea factible por 
métodos matemáticos, físicos y computacionales. El 
sólo hecho de capturar una escena implica su simpl- 
ficación puesto que por ningún motivo estamos 
midiendo todas las propiedades de la misma. Ade- 
más de esto, al sistema de captura o Instrumento 
sensor, se le puede ver como un operador definido 
por su función de respuesta al impulso. De tal forma 
¡entonces que el instrumento es un transformador de 
información que simplicala escena yla presenta en 
un formato adecuado para su manejo, interpretación 
y modelado fsico ulterior, En su expresión más gene- 
ral, esta transformación consiste en el mapeo (fig 
52,) de un conjunto de funciones de entrada (y), 1 
= 1.2,..n, aun conjunto de funciones desalidag,xy), 
|= 1,2...m, definidas en los dominios D, y D respec 
tivamente y, que se suponen ser bien comportadas 
desde el punto de vista del cálculo diferencial e 
integral. Sea ahora un conjunto de operadores 
establecidos de acuerdo al instrumento y al experi- 
'mento que se realiza para hacertal mapeo, Entonces 
las funciones de salida se relacionan con las de 
entrada de la siguiente manera 
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Una combinación lineal, una convolución y una 
Composición mukiespectralde funciones, son algu- 
nos ejemplos representativos del proceso de 
creación de imágenes dado por (5.35). Debemos 
enfatizar aquí quelos operadores “, nosólo pueden 
caracterizarlatransformación escena _ imagen, sino 
también la transformación imagen - imagen. La 
modelación matemática delos operadores ", puede 
tomarmuchas formas asícomo su medicióno estima- 
ción física; vamos a emprender aquí el camino de su 
modelación matemática sin dejar deladolas conside- 
raciones físicas pertinentes, Una manera de estable- 
cerun modelo para ”, esa partir de operadores de 
singularidad dada la connotación física que éstos 
tienen. De entre todos estos operadores, tomemos 
por su simpleza y facilidad de manejo al operador 
deltade Dirac, junto conla interpretación que tiene de 
seruna gaussiana infinitamente estrecha. Además, la 
delta de Diractiene mucho que ver con el proceso de 
muestreo de una escena, ya que éste puede verse 
como una aplicación sucesiva del operadorinstru- 
mento a través de un método de medida. Recorde- 
mos la importancia que tiene todo esto puesto que el 
muestreo de una escena es el principio de su análisis. 

Entonces, comencemos por revisar algunas de las 
propiedades de la delta de Dirac que están relaciona- 
das con el presente análisis. La función deta de Dirac 
es un operador que en dos dimensiones se define 
como sigue 


= cuando x 





0 de otra mas 


Existenotras definiciones más formales” que ésta, 
paco con la relación (5.36) es suficiente para nuestros 
propósitos. Consideremos, por tanto, las siguientes 
propiedades fundamentales 


[| scriasay 1. v0>0 


donde < es unlímit de integración infinitamente 
pequeño pero mayor que cero. Lo que estamos 
expresando con (5.37) es precisamente a una 
¡gaussiana tridimensional infinitamente estrecha pero 
cuyo volúmen bajo su superficie es igual a 1. La 
siguiente relación integral se conoce como la propie- 
dad de Shifting y se escribe como 


[ eccantéce rondgón = fa) 
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Esta igualdad es básica para entender el formalis- 
mo del proceso de captura de una escena, puesto 
¡que detodos los valores paralos cuales está definida 
lafunción, se selecciona uno de ellos: el que corres- 
ponde precisamente al punto (x,y) y, esto constituye 
la esencia de la captura de información bigimensio- 
nal. Otra manera de ver esta relación (5.38), es decir, 
que una función de entrada al sistema de imágenes 
(535) puede ser representada por una suma de 
funciones delta con factores de pesorí-Sielsistema 
de coordenadas al que se reñere la función delta es 
ú'ortonormal y si no hay correlación observacional de 
uneje de observación con respecto al otro, entonces 
se tiene 


sm 


donde las funciones deta uni-dimensionales satisfa- 
cenlas versiones en una dimensión de las propieda- 
des arriba mencionadas. De esta última relación es 
claro que cualquier propiedad que se cumpla en una 
dimensión también se satisfará en dos dimensiones. 
Un caso particular del sistema de imágenes (5.35) 
arriba mencionado y de especial interés porque per- 
mite visualizar el proceso de captura de una escena, 
es aquel donde setiene una relación unívoca entre las 
funciones de entrada y las de salida 


2) AL) 
5.40) 


donde además se cumple el principio de superpo- 
sición, es decir el operador_ es lineal aditivo. Poca 
generalidad estamos perdiendo al tomar de aquí en 
adelantelarelación dadapor(5.40), ya quela mayoría 
de los sistemas que se presentan en la vida real son 
Iinealesadíitivos, obien los diseños experimentales se 
realizan para que cumplan contal propiedad. Cuando 
no se puede tener ni lo uno ni lo otro, entonces es 
factible en una buena cantidad de casos, la aplicación 
de técnicas de quasiinealización”, Esta Última rela- 
ción (5.40) en combinación con la propiedad de 
Sifting, produce la siguiente caracterización matemá- 
tica para una imagen continua 


an 0) | empate macia 


ron) 


Esta expresión cuantifica el mapeo representado 
gráficamente por la figura5.2. La ecuación (5.41) nos 
ice que seleccionamos un punto de la escena f por 
medio dela función  ylo enviamos al espacio de 


la imagen q. Desde luego la función de entrada 
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detransformación / pueden ser invertidos en elorden 
de aplicación, por lo que se tiene 








lo que quiere decir que la imagen g es formada porla 
suma delas contribuciones de todos los valores de la 
función ftransformada por 1 . En forma ideal, cada 
punto de f debería de contribuir a sólo uno de g, pero 
'en un sistema real cada punto de f y sus respectivos 
vecinos contribuyen a formar los puntos de q, de tal 
manera que el proceso de formación de una imagen 
¡es un proceso de vecindad y no uno puntual. El que 
estosea asíse debeno sólamente ala respuesta finita 
del sistema de captura, sino también a los mecanis- 
mos de degradación que intervienen inevitablemente 
en el proceso de formación de una imagen. A estos 
dos aspectos se suma el hecho del tamaño finito del 
CIV que, como vimos anteriormente, unefec- 
10 de fitro pasa bajas. Todo esto trae como conse- 
cuencia que la imagen pierda con respecto a la 
escena: 1). contraste, 2). definición y 3). agudeza; es 
decir la imagen puede aparecer como un efecto de 
fuera de foco con poca tonalidad. Es necesario decir, 
sin embargo, que muchos sistemas modernosson de 
gran eficiencia y que portanto laimagen resultante es 
de gran calidad ya que los tres efectos mencionados 
¡son mínimos. Por otro lado es fácil ver que el equivar 
lente del modelo dado porla(5.41) paraslcaso digital 
es 
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¿con kLm y n enteros y (242.5) las coordenadas 
del centro del CIV. A partir de (5.41) y (5.43) es claro 
que podemos obtener modelos finales para 
imágenes continuas y discretas. 


Ahora bien, si el operador 2 depende sólo de las 
coordenadas (xy), es entonces un operador lineal 
Como se supuso anteriormente y entonces la imagen 
9 se escribe como 











La combinación de operadores 13) describe pre- 
cisamente la forma como se realiza el mapeo del 
espacio dela escena (1g.5.2)alespaciode laimagen, 
ss por esto que a 


Em) 0 





se le conoce como función de transferencia cuando 
se habla en términos generales de la transformación 
de una imagen a otra; sean continuas o discretas. En 
slcaso particular de unsistema de captura através de 
un sistema sensor, ah se le conoce como función 
respuesta al impulso o función de dispersión 
punto. Veamos fisicamente el porqué de estas útimas 
denominaciones para h. Con anterioridad se anotó. 
que el proceso de formación de una imagen es de 
vecindad y no puntual. Imaginómos que logramos 
enviar, al sistema de captura caracterizado porh, un 
trende pulsos deluz cuyaintensidades ly muy cortos 
en duración, tal que podamos considerarlos infinita: 
mente estrechos. Siel detector fuese perfecto, produ- 
cia ala salida un tren de pulsos de voltaje también 
infinitamente estrechos. Como nunca es este el caso, 
el detector genera pulsos con una cierta anchura y 
cuya equivalencia en intensidad es menor la de los 
pulsos de entrada. En otras palabras el tren de entra- 
da ha sido degradado por la eficiencia finita cel 
detector. 








En la fig. 5.8 se esquematiza en una dimensión 
dicho proceso de degradación y cabe anotar que 
éste, generalmentese modelacon buenos resultados 
¡como una gaussiana cuya varianza es proporcional a 
larespuesta del detector, pensemos ahora que cada 
Punto (xy de la escena es un pulso cuya intensidad 
es flx.y), ala salida de h, este punto estará represen- 
tado por una gaussiana tí-dimensional. Para un 
punto y su vecindad, las respectivas gaussianas se 
traslaparán más o menos dependiendo de h, por lo 
'queelvalorde un punto delaimagen será el esultado 
dela interacción (a través de h) de una colección de 
puntos de la escena. Entre más anchas sean las 
'gaussianas más puntos de la escena contribuirán a 
uno de la imagen. En principio las colas de las 
gaussianas son de infinita extensión, detal forma que 
todos los puntos de la escena contribuyen estricta- 
'mente 2 cada uno de la imagen; en realidad esto no 
esasícdadala escalafinita de valores (2") que seutiza 
para cuantificar los valores de los puntos de la ima- 
gen, porlo que más alláde un ciertolmitelas colas no 
contribuyen a la formación de la imagen. Los meca 
nismos de degradación? pueden ser explicados 
tambiéncomo unas de gaussianas yde 
aquí que cuando hay este tipo de interferencia en la 
¡generación de una imagen, podemos pensar que la 
contribución a un punto de ésta provenga de una 
vecindad mayor. El razonamiento anterior es clara- 
mente aplicable a la formación de imágenes conti- 
nuas para cuando el sistema de captura es Óptico; 
¿peroqué pasacon imágenes digitales? Alhacerse la 
integración sobre un elemento CV, es claro que un 
pixel de la imagen digital proviene de un conjunto de 
Puntos de una vecindad determinada por tal slemen- 
to, ya que éste es un subconjunto compacto y conve- 
xo de la escena. Los mecanismos de degradación, 
porotrolado, hacenque CIV vecinos contribuyan ala 

'de aquel quese encuentra bajo medición. 
Todo lo anterior nos indica que en cualquier caso se 
tiene un proceso de vecindad además de haber un 
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Fig. 6:8.- Tren de pulsos de luz muy cortos en duración y respuesta del detector h representada por el 
tren do gaussianas a/a salida de h. Por esta respuesta finita, la intensidad de los pulsos se ve disminuida 


en una cierta Al. 
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mecanismo decorrelación observacional, es decirlos 
elementos de laimagen digital están correlacionados 
por el método de medida, inelucibiemente impuesto 
por el sistema de captura, Desde luego aparte de 
sto, siempre hay una correlación entre valores de 
pixels vecinos debido a la variación natural de los 
valores de radiación que provienen dela escena. Es 
la correlación observacional la que eventualmente 
hayque eliminar antes del análisis delaimagen digital 
correspondiente, 


Volvamos ahora al desarrollo del modelo para una 
Imagon; la combinación de las expresiones (5.44) y 
(5.45) nos da el siguiente resultado 





para una imagen continua. Para cuando el operador 
no es lineal, se tiene entonces 


es decir, la función de transferencia depende de los 
valores de la escena. Un ejemplo fípico de esto es la 
respuesta no-ineal de un fim fotográfico en los extre- 
'mos de su rango dinámico. En sistemas lineales, la 
descripción de la función h en términos de las cuatro 
coordenadas (x,,.:.,) esla más general posible y 
le permite variar con respecto ala posición tanto en el 
plano objeto (escena) como en el plano imagen. Sin 
embargo, muchos sistemas de formación de 
imágenes digitales, actúan uniformemente a través 
de los planos mencionados de tal manera que para 








ste caso lineal, se tiene una expresión de convolu- 
ción entre hy 





zar fisicamente esta convolución imaginando que 
sistema de captura "analiza" en un proceso de barí- 
do, elemento a elemento de la escena a través del 
plano '“.1) ; veamos esquemáticamente esto por 
medio de la sucesión de figuras 5.9a-e. En la primera 
de esta sucesión de figuras se aprecia la extensión 
fisicade una escena cualquiera 1%.) y enla segunda 
sohadibujadouna función detransferencia 4.1) . 
Aparticde aqui sigamos larepresentación geométrica 
de una convolución; así en las figuras 5.96 y 5.9d se 
han dibujado hí6.-0) Y hix-t7-1) respecti 
mente, para finalmente concluir en la figura 5. 
donde se muestra en achurado, el traslapo de h y 1 
Hay que recordar que tanto h como f son superticios 
yque enlas figuras 5, 9a-e sólo se ve un corte de éstas. 
Lazonasombreadade la igura5, Ses el producto de 
hy fque a su vez es una superficie. El volumen bajo 
esta superficie es precisamente la doble integral de h 
yt. enotras palabras, este volumen representa elvalor 
dela convolución en (x,y). Como este último valor es 
el que se le asigna a la imagen en dicho punto, 
podemos concluir lo siguiente, recordando siempre 
el modelo digital dado por (5.43): El sensor remoto 
del sistema de captura subtiende un CIV cuya exten- 
sión física es igual a . La manera como este 
sensor responde a la radiación que proviene de él, 
está dadapor  1it.1:=22 0) , donde (mun) son las 
¡coordenadas de su centro, El volumen mencionado 
representa por tanto la cantidad de radiación 
proveniente de la escena y subtendida cada vez por 














el sensor sobre una subárea (CIV) de ésta. Al 
discretizar el tamaño de este volúmen tendremos el 
valor del pixel corresponciiente, de la forma como fue 
explicado en la sección 5.2. Con las consideraciones 
arriba planteadas es muy fáci escribirla versión digital 
del modelo matemático de creación de una imagen 





con (M,N) la dimensión de una imagen digital 
rectangular. 


En las relaciones (5.48) y (5.49) se han plasmado 


las transformaciones imagen continua — imagen 
continuas imagencontinua — imagendiscreta. Para 
fines prácticos Í( £(L.n) — ) puede ser vista como 


una imagen continua y de esta forma tenemos el 
'mismotipo de modelo para la capturade una escena, 
¡como para latransformación de una imagen. Esto es 
posible porque no se ha hecho ninguna suposición 
explícta que diferencie una escena de una imagen 
continua. Nos falta escribir todavía la transformación 
imagen discreta. = imagen discreta. Al suponer que 
1 está definidasólo para valores discretos de 





tenemos rápidamente el modelo requerido 


elit) EY Ennio 





donde hacemos notarque — E - 
ya que my nson discretos y portanto 2 = 
, ambos iguales a uno, Desde luego cuando 
hayinvarianzafrentea translaciones, (5.50) se escribe 











ua 


'Cabeagregar ahora quelarelación paralatransfor- 
mación imagen discreta - imagen continua, sólo 
puede darse a la luz del Teorema de Muestreo. 


Enunaespecie deresumen consideremos ahorala 
tabla5 2 y lafigura 5.10, en donde podemos apreciar 
esquemáticamente los modelos y consideraciones 
úhastaahora desarrollados. La figuraS.10 está relacio- 
"nada con la 5.2 pero describe en forma más concisa 
la transformación de una imagen. En la tabla 5.2, se 
han escrito los principios de causalidad y conser- 
vación de energía para cada uno de los tres modelos 























Fig. 5.9.- Representación esquemática de la conwolución de h y f para producir la imagen 9. 
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Fig. 5.10.- Representación simbólica de los diversos modelos de transformación de imágenes. 




















establecidos. El primer principio nos dice que no hay 
energía radiante negativa y el segundo puede verse 
'comounaleyde conservación de lainformación. Esto 
último siempre y cuando se cumplan las reglas del 
teorema de muestreo bajo condiciones reales de 
operación del sistema de captura. Nos falta agregar 
una cuestión de nomenclatura: diremos que cuando 
estamos capturando a una escena, el instrumento 
correspondiente se caracteriza pora función de res- 
puestaal impulso, cuando setransforma en general a 
tuna imagen ah se le conoce como función de trans- 
ferencia, pero si la transformación se ejecuta en el 
dominio espacial entonces a he lellama el Kernel de 
la translormación. De esta manera estamos nom- 
brando ahdetres formas diferentes dependiendo del 
contexto bajo el cual se transforma la imagen, aún 
cuando se utilice, como ya se describió arriba, el 
mismotipo de modelo matemático detransformación 
deimágenes. Es necesario dejar claro que los mode- 
los planteados son aplicables para cuando se genera 
una imagen muhiespectral o bien para cuando se 
Opera sobre ésta. Baste recordar que se había intro- 
'ducido enlasección 5.1 un parámetro para connotar 
intervalo de energía, entonces una imagen digital 
Muhiespectral estará descrita por 
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donde : esfinalmente un conjunto de intervalos de 
energías, no necesariamente contiguos ni traslapa- 
dos. Á cada imagen para un intervalo, se le conoce 
como banda y al conjunto de imágenes para un 
conjuntodeintervalos, sele llamaimagen multiespec- 
tral. 


Elmuestreo de una imagentrae consigo lacaracte- 
rización matemática de la misma en forma digit 
Puesto que el muestreo de información necesar 





mente implica su discretización, es decir que al 
muestrear ordenadamente una escena se produce 





solamente ía reconstrucción adecuada de la infor- 
mación, sino también para minimizar la creación de 
artefactos inherentes a todo proceso de medida. La 
reconstrucción se refiere a las técnicas de ftrajo 
adecuadas para eliminarlos artefactos que de todas 
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maneras se introduzcan y ala interpolación necesaria 
para obtener una imagen continua que represente a 
la escena con la mayor fidelidad posible, Luego 
entonces el proceso demuestreo dela escena genera 
una función discreta que es necesario estudiar con 
detallo, Para esto diremos que una función discreta 
estádefinida solamente para valores discretos de sus 
variables, por lo que la representaremos como una 
secuencia de números enteros, Además de esto, la 
función misma sólo puede asumir valores enteros 
también, El resultado es un conjunto ordenado de 
números que representan las muestras obtenidas, 
tanto espacial como radiométricamente, de una fun: 
ción analógica o continua. Cómo se obtienen estas 
muestras y cual es su representación matemática es 
justamente lo que se ha discutido hasta ahora. Dire- 
mos portanto que una función discreta es un conjun- 
10 ordenado de números que representan los valores 
deuna escena para lugares (IV bien especificos de 
a 

















dondelos corchetes (| representan una secuen- 
cia ordenada de números enteros. Aquí recordemos 
quela discretización espacial sobre (xy) y que produ- 
¡enteros asociados alos pixel (na es 
lo llama igitización. A la discretización r 
sobre g y que produco la mencionada secuencia sele 
llama cuantización. El conjunto de intervalos de ener- 
gía (os 1, está definido por el diseño 
Particular del sensor mutiespectral. Los números 
(n,.n,) serán considerados también como números 
naturalos, poro pueden ser vistos como enteros sl el 
rigen del sistema de coordenadas se fja:en el centro 
de la imagen. En consecuencia los números (,.n¿) 
pertenecen al siguiente conjunto 





















donde Y enelgnipo dejos enteros yla imagen 
os cuadrada N' lado el proceso de 
diosatcación de ls “rioraiada ergía que 
provienen delaescena, requiere que cadanivel delos 
Pixels ostó restringido a números enteros de m-bits y 
esto desde luego con propósitos de manipulación 
úficionteencomputadora de laimagen digital. Asíque 
finalmente los valores de los pixels de cada banda 
¡son miembros del siguiente conjunto 
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donde | es el sem-grupo de los números natural 
Los números (n,, n,X) pueden desde luego tomar 
sólovalores discretos y enla prácticam = 7.u8,lo que 


significa que el rango para los niveles de los pixels de 
laimagen muliespectral, estárestringido entre 0 y 127 
Obienentre0y256 respectivamente. Estees un rango 
conveniente porque cada pixel ocupa un byte de una 
computadora digital. Hemos representado hasta 
ahora los elementos (pixels, bandas) de una imagen 
digital multiespectralutlizancio modelos determinísti- 
cos yhemos dejado deladoel aspecto estocástico de 
la misma porque aún no la hemos visto en su conjun- 
to, Enlasiguiente sección estudiaremos la naturaleza 





estadística de la generación y transformación de 
imágenes digitales. 
5:5- Caracterización Estocástica de una 





Imagen Multiespectral 


Para distinguir claramente entre los elementos 
deterministas y los estocásticos de una imagen digi- 
1al, os conveniente establecer ahora una clasificación 
de cuatro clases de imágenes, Esta clasificación 
cualitativa tiene que ver con la representación visual 
de una imagen pero también ostá relacionada con la 
manera on la que se le procesa y se la representa, 











división entre clases no está perfectamente definida 
sino más bien existo ciortotraslape entro llas, aún así 
la clasificación es suficientemente útil para nuestros 
propósitos. La primera clase es la más amplia y es la 
que comprende a las imágenes más compleja 






das, Aquí encontramos todas las imáger 
"nas comunes en la vida real: la foto ce un palsaj 
retrato de una persona, una radiografía y una imagen 








ocurren normalmente en la naturaleza y no de esce- 
nas anificialmente creadas por el hombre, En dond 
interviene la generación de imágenes por algún 
medio mecánico, tipográfico o electrónico, con es- 
tructuras relativamente simples y pocos tonos, es en 
jundaclase. Aquíse situan los 
gráficas y las pinturas, representados en 
forma digital. En esto grupo de imágenes las clases 
de cobertura están bien dolineadas y representadk 
por un sólo color. La frontera de las estructuras 
comprendidas está normalmente bien definida. 
¡Como podrá apreciarse la diferencia en complejidad 
ntro la primera y la segunda categoría de imáganas 
ds bastante grando, Es raro encontrar un caso que se 
"encuentre en la frontera de estas categorías. Algunas 
escenas naturales muy simples, como un patrón de 
aperiónci nextrones, podian tene una repasan. 
tación digital que talvez cayesen en la segun 
deimágenes. Sin embargosi conocemos el origen de 
la imagen y aplicamos la regla de selección de la 
primera clase en el sentido de que la imagen no debe 
provenir de medios aríficialos, entonces desaparece 
la ambigúedad. Sino conocemos el origen de laima- 
¡gen en cuestión, la confusión puede prevalecer y una 
decisión arbitraria será necesaria. La tercera catego- 
ría comprendo imágenes más sencillas aún, pero no 
es muy diferente de la segunda categoría. Aquí se 
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encuentran aquellas imágenes que están compues- 
tas porcontornos solamente formados porlíneas que 
pueden serde varios colores, Es comosi de un mapa 
hubiésemos dejado sólo la frontera entre regiones. 
Deestaformasi enla clase dos ponemos un mapa de 
los países de América, en la clase tres podría quedar 
la imagen de los contoros de las fronteras entre 
dichos países, cada frontera con un color diferente. 
En esta tercera categoría podemos también incluir 
cieno tipo de gráficas sencilas, En la cuarta y última 
categoría se encuentran todas aquellas imágenes 
binarias, compuestas por puntos o polígonos de un 
¿olor y un fondo uniforme de otro color. Un ejemplo 
'muysocorrido de estaciase deimagenes es dondese 
tienen contomos en blanco y negro únicamente. La 
imagen de una página de un libro consttuye un 
ejemplo adecuado de esta última ciase. Aquí caenlas 
imágenes estructuralmente más sencilas que pue- 
den ocur. De la clasificación ariba establecida 
podemos concluir varias cosas. En primer lugar 
podemos afirmar que las gráficas forman una clase 
de imágenes y que pueden ser operadas como tales, 
si es necesario, con la misma herramienta matemáti 
ca que la empleada para las imágenes más compie- 
jas. Así, por ejemplo, la gráfica de la distribución de 
intensidades del campo magnético terrestre, puede 
serfítrada y realzada con consecuencias muy intere- 
santes desde el punto de vista geofsico, En segundo 
lugary hablando entérminos generales las imágenes 
clase 2, 3 y 4 pueden ser descritas por modelos 
deterministas debido a su relativa sencillez ya que en 
muchos casos son propiamente generadas por 
modelos matemáticos preestablecidos. En tercer 
lugar,lasimágenes querequieren de una descripción 

ica son aquellas que se encuentran en la 
clase uno, siendo las más complejas ylas que apare- 
cen normalmente en la investigación experimental. 
Son estas imágenes las que por lo común son de 
carácter multiespectral. En cuarto y Útimo lugar, es 
claroquealoperar una imagen de entrada, através de 
las relaciones 5.46 6 5.49, se obtendrá una salida de 
la misma clase que la de entrada, siendo además la 
descripción respectiva del mismo!ipo: determinista 
estocástica. Esto se debe a que el operador h es de 
tipo determinista. Entonces las imágenes de ciases2, 
3y 4norequieren de unamayor consideración detipo 
estadístico, además de que existen una buena varie- 
dad de tratados sobre gráficas que se avocan espe- 
cíficamente a esta clase de imágenes. La clase uno sí 
tequiere de consideraciones estadísticas y para esto 
es necesario hacer una revisión de algunos 
conceptos relacionados con variables aleatorias y 
procesos estocásticos, Como no es el objeto de este 
trabajo el estudiar de por sí estos problemas, la 
revisión que se hará será en el contexto de las imáge- 
nes digitales muliespectrales. 

Considérese ahora un experimento de naturaleza 
sstadística*"%, es decir, uno donde las respectivas 
salidas están regidas por leyes probabiísticas. Como 
un ejemplo de esto pensemos en seleccionar al azar 
una fotografía de entre un conjunto muy grande de 
imágenes. Esto teniendo en mente que para imáge- 
nes de digamos 1024 x 1024 elementos, podemos 


eventualmente reproducir casi cualquier escenade la 
vida diaria, si contamos con un número suficiente: 
mente grande de pequeños cuadros de 128 posibles 
tonos de gris. Para imágenes en color el número es 
tres veces más grande, pero por el momento even- 
tualmente reproducir casi cualquier escena de lavida 
diaria, si contamos con un número suficientemente 
¡grande de pequeños cuadros de 128 posibles tonos 
de gris. Para imágenes en color el número es tres 
veces más grande, pero por el momento el ejemplo 
contonos de grises adecuado. Es mucho más di 
entender experimentos de esta naturaleza para imá- 
¡genes que para series detiempo, puesto que setiene 
una interpretación visual de las respectivas salidas y 
porque éstas son muchisimo más numerosas que en 
una dimensión, Esclaramente más fáci entender que 
alrodarun dado estamos ejecutando un experimento 
denaturaleza estadística y que al observar una de las 
caras del mismo, estamos en realidad observando 
una de las seis posibles salidas que existen con igual 
probabilidad cada una. En dos dimensiones debe- 
mos considerara laimagenen su conjunto, comouna 
posiblidad de entre muchas que existen: 
¡2s%Hsas donde P significa permutaciones. 

Sea ahora * el conjunto de posibles salidas del 
mencionado experimento, que como ya se explicó 
arriba, puede llegar a ser muy grande 








Paraque < esté bien definido se requiere que 
las condiciones delexperimento estén determinadas, 
es decir que ss encuentre fijada la goomatría sensor 
Objeto-fuente de radiación y que los rangos de medi- 
da de las variables involucradas hayan sido conve- 
nientemente dados. Unasalida -, consiste portanto 
en la selección de una imagen de entre un conjunto 
muy grande de ellas. Además de esto asignamos una 
cierta probabilidad p, a cada salida =, < 2 ,detal 
forma que 





loque quiere decir queno puede haber unasalidacon 
probabilidad negativa, ya que físicamente por al 
Principio de causalidad no existen experimentos con 
E, +8. Consideremos un subconjunto / de 
talque 4:  ,entonceslas salidas asociadas a 
tienen también una probabilidad bien definida, A 
los subconjuntos , de £ selesconoce como 
eventos y enla práctica es necesario considerar estos 
eventos y sus probabilidades asociadas. Puede 
ocurrir en la realización de un experimento que un 
evento sea igual a una salida aún a“; algunos 
ejemplos de esto son 
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Estas últimas consideraciones nos dan pauta para 
ver la cuestión de la captura de una imagen de la 
siguiente manera: podemos considerar las salidas 
««, con cierta fexibilidad en el momento de discret 
zaruna escena, de tal forma que la generación de un 
Pixel a la hora de medir la radiación sobre un CIV, 
Puede ser considerada como la salida de un experi- 
mento que genera pixels para conformar una imagen. 
Así, entonces, una imagen estará constituida por un 
conjunto ordenado de salidas (pixels). Viendo las 
cosas de este modo, 1 está formada por tanto por 
todos los pixels posibles de una imagen. De aquí 
podemos afirmar que, ll, 0s,--%, representa 
el conjunto total de imágonos quo pueden sor 
conformadas para una geometría y condiciones 
específicas de un experimento de naturaleza 
estadística. Una imagen multiespectral producto de 
este experimento, estará formada por un 
subconjunto de 7... Poratro lado cuando se barre 
una escena con un barredor multiespectral para 
obtener una imagen digital, se obtiene un cierto 
conjunto de pixels que forman cada línea de barrido, 
en esto se produce un subconjunto de 5: A = 
=== %a!, con lo que se ve claramente que una 
línea de barrido puede ser vista como un evento 4; 

“con tantas salidas como pixels tenga la dimensión 
dela imagen en la dirección de barrido. Habrá tantos 
eventos 4, como líneas tenga laimagen. La imagen 

, estará formada por el conjunto %, + (Mis, 

¡que constituyen el número total de lineas de 
barrido de ósta. Por ejemplo una imagen 
multiespectral de — bandas y dimensión espacial 1- 
líneas x m-pirels, será producto de m eventos con 1 
salidas cada uno y cada plano de una banda estará 
representado poruna 4, porloquetendiemos 

y deestas  , - Esnecesario puntualizar aquí que 
de acuerdo a lo dicho en este contexto se tendrá 

Y = k, ya que para una geometría dada podemos 
¡generar muchas imágenes muliespectrales de una o 
varias escenas siempre y cuando éstas se ajusten a 
las condiciones del experimento. Entonces, en 
general k podrá ser muy grande y + asumirá valores 
relativamente pequeños yaque + serefiere a unasóla 
imagen multiespectral de las muchas que pueden ser 
producidas. Con esto hemos generalizado lasefin; 
iones arriba dadas (5.56 y 5.57) ampliando nuestro 
crterio sín contradecir las consideraciones estadíe 
cas originales. 

Ahora podemos agregar que un evento puede 
ocurrir cuando cualesquiera de las salidas que lo 
componen ocurra. Un evento ciertoes aquel que está 
compuesto por todas las salidas posibles que com- 
ponen una banda de una imagen multiespectral, así 
por ejemplo 1; es un evento cierto. Los eventos 
cumplen con propiedades de la teoría de conjuntos, 
así por ejemplo/, es el complemento de — 4, 

y constituye la no-ocurrencia de 4, , además 

A, 18,7 4,0 A, son también 
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eventos. A 1. sele conoce como el evento nulo y es 
por defiición un conjunto vacio,  entoncessi 
y son mutuamente 
pelle pa Jeremos ahora la probabili- 
dad P4) asociada al evento / , detal manera que 
setiene que 





Ahora bien si 1 , es una secuencia de 
eventos mutuamente excluyentes, entonces 








esta secuencia bien podría ser la sucesión de líneas 
de barrido en el proceso de creación de una imagen 
digital. Por otro lado, si a cada salida “1 de un 
experimento asignamos un número £(%,) —,lo que 
es entonces una función sobre el conjunto. %, = 
Lo =.! y que recibe el nombre de variable 
aleatoria, La relación fentre las salidas “y 
ylos números — (4) ds doiipo determi: 
rista, ss decir, unavez que ha ocurrido una salida “y 
el valor quetoma la variable aleatoria es una cantidad 
completamente determinada a través de la regla de 
asignaciónf — .Loqueesaleatorio eslaocurrencia 
dela salida misma, la que no es un número sino un 
hecho que ocurre, Es muy importante señalar que de 
acuerdo a lo discutido hasta ahora, la regía de 
asignaciónt « estádeterminada por elsensorremoto, 
, constituye el hecho de medir sobre un CIV 
inado y *(u;) es el número que se le asigna 
| correspondiente. Como yase había apuntado 
en lasección (5.2), enel momento de barrer una línea 
sobre la escena, se produce una serie de tiempo de 
los valores de radiancia provenientes de la escena; 
dichaserie constituye portanto una variable aleatoria. 
Existen otros tipos de eventos que son interesantes 
paralos propósitos de andlisis de imágenes que nos 
hemos fado. Asi, por ejemplo, el evento. 11 3 «1 
representa el conjunto detodas las salidas “, para 
las cuales — 1(4/) <<, dondezesunnúmero 
arbitrario pero fo, de tal forma que para cadaz, 
(Les un evento, También, para dos 

















números 21 yea don 2) esunevento 
formadoportodaslas salidas paralas cuales “(4,) se 
encuentra entre +, yz» .Parafinalizar en esta 


serie de ejemplos, es claro que el conjunto (+= A 
es un evento en forma similar a los demás. De estos 
eventos tomemos por conveniencia el conjunto 
1 + +1, la razón de elegir a éste se verá en el 
desarrollo subsecuente. La probabilidad asociada a 
es un número que depende de 2, es decir 
es una función de z. Sea esta función P.(2), entonces 


Po 
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Ala función P/(z) se la conoce como la función de 
distribución de probabilidades de lavariable aleatoria 
1. Supóngase que Pz) es derlvable y que tiene las 
siguientes propiedades 





a) MO Pt estaúlimacsel 
¡conjunto de todas las salidas yla pimera so debe al 


principio de causalidad. 


Dj. Pas) 5 PCs) paro 21 5 7, ¡esto.esasíya 
que para 24d , laa) es 
siempre un subconjunto de — (£ <=: yademás 


si ACB, entonces Ps 10 


O): Mes 8 18 29) = Be(2i)= F(a1) que se deriva fé: 
cllmente delo anterior. 


Podemos ahora introducir la siguiento cofinición; a 
la dorivada 
0) 
rm 


(5.62) 


de la función de distribución se la llama función de 
densidad de probabilidad de la variable aleatoria. 
Con la combinación de esta únima definición y de la 
propiedad (c) puede obtenerse un resultado útil con 
luna interpretación sencilla 


Mas al Pa 


Po 


il 
pta 
Y 





0.6) 


con lo cual es claro que 


(50) 


Con este último resultado podemos entender fácik- 
mente que la definición (5.61) es adecuada, pues se 
trata de una probabilidad acumulativa ya que se 
consideran las salidas para las cualos (4) E 
al aumentar este número 2, se incluyen todas las 
salidas anteriores más las que se agregan on esto 
incremento, En otras palabras y ala luz de (5.64) so 
considera que el área bajo la curva de la función de 
densidad es la probabilidad asociada al evento 
(£ s 2) .A partir de esta interpretación puede verse 
que cuando el conjunto de muestras disponibles es 
representativo del universo completo de medidas, 











entonces setienequesi 1 essuficientemento 
, se obtiene la siguiente aproximación 
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6.65) 
de aquí que la función de densidad pueda escribirse 
como 


PC) La MR ea d)/e 
Ao 


Estadefnición tiene un valor práctico yaque permi- 
o visualizar un método aproximado para obtener la 
función de densidad tanto de una variablo aleatoria 
continua como discreta, Este método es el siguionto 


1.- Ropitase el experimento un gran número N de 
voces. 


2. Cuéntese el número de veces que la variable 
Aleatoria asumo valores entre » y + +4: ,conz 
« << al intervalo de valores que toma la variable 
aleatoria; sea este número 1, 


3, Es claro que la probabilidad de que la variablo 
aleatoria tome valores entro < y ; + 32 es de apro- 
ximadamente ', A 

4. Portanto», 





Estoresuttado refuerza la conveniencia do ladofini- 
ción adoptada para P(), ya que del punto (4) e dico 
intoncas que, para una banda de una imagen digital 
muliespectral, p(2) está representada, on buer 
aproximación, por el histograma de valores de los 
pixels que la componen. Nótese que el área bajo la 
Curva de la función de densidad entra z, y z, os igual 
ala probabilidad de que ftome un valor entro 2, y 2,. 
Do aquí podemos decir que +1») > 11002) 'esía 
probabilidad de que ftome un valor entre > y +, La 
valuación, modelacióno estimación dela función de 
densidad nos permite hacer una descripción de 
Imagen por sus características promedio, es decir, de 
la imagen en su conjunto. Dos de las propiedades 
promedio más interesantes son el valor medio y la 
desviación standard. El valor medio o el valor 
esperado E HI£) > f de una variablo aleatoria so 
define como 
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(5.07) 
Bajo ciertas condiciones!? , , que se conoce 

también como el primer momento estadístico, 

igual o aproximadamente igual al promedio 





A 
e (5.58) 
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entonces utilizando las anteriores definiciones, se 
tiene 
HO Ta Mr, 3 5) 


expresión que puede ser reescrita como sigue 





que es la definición usual de promedio. En la relación 
(6.70) se supone implícitamente que se cuenta con 
una muestra válida del universo de medidas. La 
¡segunda propiedad promedio utlizada en la descrip- 
ción del conjunto de una imagen, es la desviación 
standard o segundo momento, Una manera intere- 
santo de introducir esta cantidad es a parir de la 
misma relación (5.57) que se empleó para el valor 
medio, puesto que es fácil demostrar que si 4 = 1111 
+s una función de la variable aleatoria f, entonces y 
fartión ae LA vacia aatora y era la gun 
cel 
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sm) 
quees en realidad una generálización dela expresión 


(6:87). Sea ahora “: el valor medio de , entonces 
lavarlanza 0!  def,se define como sigue 


up tocara 
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Alaraíz cuadrada dela varianza sele conoce como 
la desviación standard y también existe una relación 
aproximada en forma análoga a como lo es para el 
valor medio. Existen otros momentos con interpreta- 
ción simple y que se emplean en la descripción de 
imágenes, sin embargo, éstos se verán mejor en los 
siguientes párrafos donde se generalizarán las ideas 
vertidas en la presente sección. Los momentos de 
alto orden no se tocarán aquí pues están fuera de los 
objetivos de estaobra. Para lavarianza también existe 
unainterpretación sencilla como lo es que representa 
na medida de la dispersión de los datos o una 
medidadel contenido delaintormación. Alproporcio- 
nar la función de densidad de una imagen estamos 
finalmente dando una descripción de conjunto de la 
misma, yaquenonadamás damos la probablidad de 


ocurrencia de un cierto valor de un pixel, sino que 
también a partir de dicha función de densidad pode- 
'mosobtenerlos momentos mencionados, La función 
dedensidad es portanto un modelo estocástico dela 

imagen. Existe una gran variedad de formas para le 
fncón de densidad, algunas de ells se lstan a 
continuación 





La densidad Gaussiana 
ID) = roja) exp LO al, 


61m 





La densidad de Rayielgh 


0.) 


La densidad Log-Normal 
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525) 
La densidad Exponencial 
pon sepia le, 191) 


(5.76) 
La densidad Uniforme 
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(5.1) 


En las expresiones para estas funciones de densi- 
dad, (1) es el valor del pixel cuya posición es (K,1) 
yx es una constante. Para imágenes reales es 
necesario evaluar la función de densidad para a > 0 

ya que no existen valores de radiación negativos. 
Para finalizar esta parte diremos que es importante 
hacer notar que en la realización completa de un ex- 
perímento, se tienen los siguientes aspectos del 
mismo: 1). El fenómeno bajo estudio; 2). El proceso 
de captura de información de tal fenómeno; 3). Los 
datos obtenidos; 
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Pasemos ahora a generalizar las ideas y definicio- 
nes arriba vertidas en cuanto a un experimento de 
naturaleza estadística, con el objeto de poder descr- 
Bien forma global una imagen muliespectral Sea 
por tanto una famila de funciones '.(+;' , la que 
Geponde del parámetros y quese encuenta etnia 
sobre el conjunto de todas las salidas %, - (%..4 

4,1, de un experimento dado. El parámetro s es 
talque's < 1 donde _ Les un intervalo del eje realo 
una región de un espacio multidimensional, el cual 
por simplicidad se considera Euctidiano. Esclaro que 
apanirde estas definiciones, para un valorfjo des, se 
tiene una función bien definida sobre el conjunto de 
salidas“, .Paracadavalor de *;, es decir para 
cada salida, “,(“,) es una función de s sobre un 
intervalo o unaregión del, segúnseaelcaso. Cuando 
es unidimensional y además está consituido por 
una colección de puntos discretos: : > (1.2 
entonces la familia de funciones es simplemente un 
conjunto de variables aleatorias. Existen por tanto 
dos posibilidades, dependiendo del, de visualizar a 





a).-!,(s,) es claramente una famila de variables. 
aleatorias, cada miembro de la famila generado por 
un valor de s, 


b).-'-. puede sertambién una famiña de funciones 
des, cada miembro de esta familia correspondiendo 
auna salidas, . 


¡Cuando 1 es unidimensional, la familia de funcio- 
nes se llama proceso estocástico. Cuando | es n- 
dimensional, la famila de funciones *,'-,/se llama 
campo aleatorio, el que puede ser real o complejo. 
'Supóngase ahora que | es bidimensional, es decir 
sea el plano x-y, de tal manera que s es un punto de 
ste plano ya que s < 1. Portanto s puede serrepre- 
sentado ya sea por sus coordenadas o por su vector 
deposición . Por esto, el campo aleatorio! '- será 
denotado de ahora en adelante como fr, ). de aquí 
que para un valor der. esuna: aleatoria, 
Mientras que para unasalida=, +.x-, esuna función 
sobre el plano x-y. De acuerdo a como se habían 
planteado las cosas más arriba, a *, sele había 
denominado como el conjunto de salidas correspon- 
diendo auna banda de una imagen mukiespectral, de 
tal forma que denctaremos más bien a esta imagen 
eomot(r.%), donde: es el índice que representa 
las bandas que la componen. Esta representación se 
encuentra esquematizada en la figura 5.11, donde se 
muestra una imagen multespectral de y bandas y 
caracterizada por el campo aleatorio ++. ). En este 
caso ', correspondea la selección de una imagen 
a parti de una colección de elias. En la figura 
5.11 se ve que para un valor der, esdecirr = %1, £ 
(e. ses una función de — *, ,es deciruna variable 
aleatoria. Por tanto para r fa, la variable aleatoria 
correspondiente representa los valores de un pixel 
sobre todas las bandas de la imagen y a esto se le 
onoce como respuesta especral Para una banda 
dada, *1: 1.) esuna función bi<timensional 
en el plano xy. una banda de la imagen 














Fig.5-11.-Imagenmultespectralde bandas, 
caracterizada por el campo aleatorio ((.... 
donde: = 1.2 











multespectral. Hay que hacer notar que en este 
último caso y dependiendo de como se considere a 

A, —, tíx.2,) puede representar un pixel o una 
subimagen de la banda. En resúmen diremos lo 
Siguiente: a).Para “, fia, K(r.,) es unasuperfcio 
Que representa la variación de intonsidad do la 
raciación medida enfunción de (xy) para unintervalo 
prefjado de energía ysobre el dominio D; a.osto sele 
Conoce como banda; t). Para r fo, £(z.,) es una 
variable aleatoria que se conoce como respuesta 
espectral y esto representa la variación de intensidad 
de la radiación medida para un mismo CIV y en 
función de su energía. 





¡Como ya se mencionó, el campo aleatorio tr.) 
constituye una variable aleatoria para un valor espec 
fico de 1, en general esta variable no tendrá las 
mismas propiedades estadísticas paratodos los valo- 
resder < Dparalos cuales está definida la imagen. 
Enotras palabras, lasfunciones de distribución y den- 
sidad para la famila de variables aleatorias t(r,%) 
dependerá del valor de r, por tanto, 





(5.18) 


eslafunción de distribución de lavariable aleatoria 1r 
)en el punto r. + D en el plano x-y. La función de 
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densidad respectiva está dada entonces por 





Dada la forma en la que hemos desarrollado estos 
últimos conceptos es facil verque (x., ¡está definida 
para toda '. -, poro que de ahora en adelante se 
denotaráa  <(r,%, 0 por (1) opor 1,(x1) para 
'una imagen continua o por 1! (4,1) para una imagen 
digital muliespectral. Entenderemos por tanto que 
¿ (Urepresenta una familia de funcion 
bidimensionales, Por otro lado para los puntos f,, 


Frsta;>=5, E Dl6nemos q variables aleatorias: 








..) + 4, De aquí podemos decir que una 
imagen digital de dimensión n x m está compuesta 
por q = n x m respuestas espectrales, todas ellas 
interrelacionadas a travós de la función de 
distribución conjunta 





MAA 
AS 


y atravós de la función de densidad conjunta 


pe it) 
' 








Cabe añadir que estas funciones conjuntas son 
necesarias en la descripción de la correlación espa- 
cial que existe entre diferentes objetos presentes en 
unaimagen muliespectral La modelación de (5.80) y 
de (5.81) nos da una caracterización global 

distribución de valores de radiación de la escena. 





Dado que el campo aleatorio — £(r.%,) esuna 
variable aleatoria para una r dada, entonces, en 
general, larespectiva función de densidad dependerá 
del valor de r que se escoja y por tanto el valor 
esperado de f debe también ser una función de 
entoncer 








at ee = | ensimón 





sele llama la media del campo aleatorio 





ryes asu 


vez una variable aleatoria. Un momento estadístico 
de mucho interés para el análisis de imágenes mul- 






campo aleatorio (1) está definida como el valor espe- 


rado del producto de dos variables aleatorias 1(r) y 
K(r¿) y se expresa como sigue 





TS 


Más adelante veremos que así como existen expre- 
siones aproximadas para los dos primeros momen- 
toscuandos — [esunintervalooneleje real, también 
bajociertas condiciones se pueden encontrar relacio 
nes para Ur y Ry que facilten su evaluación para 
cuando | es una región o dominio del plano x:y; Por 
lo pranto poramos decir que a mecida quer, yr, sa 
alojan entro sí, () y ((r) se hacen más distintas pues 
probablemente correspondan a objetos diferentes, 
mientras que para , y,Cercanos, es posible que f() 
y (1,) correspondan A la respuesta espectral del 
mismo objeto. Para dos objetos radicalmenta 
Mc) y (to) serán 
Y Y, Ostón relativamente 
cercanas, Una variación rápida O lenta de los valores 
delos pixels ao largo do x o y será detectada por la 
unción de autocorrelación ya que ésta es una función 
por estas razones que a Ry se lo utiliza 
como una primera aproximación para Guantificar la 
toxtura de una Imagen digital. Otro momento 
empleado para describir a una imagen os la 
autocovarianza que se define por medio de la 
siguiente relación 





























Egglrisea) = MULA) > 





Upfr ln) = Utend 


(5u0) 


La autocovarianza tiene una interpretación física 
similar a la autocorrelación ya que es muy fácil de- 
mostrar que 


Cpplen cra) e Rent) ale da (e) 


685) 


La autocovarianza y la autocorrelación están rola- 
cionadas a travós de la media “+, que bajo ciertas 
condiciones se puede interpretar como el promedio. 
De aquí que cuando estas condiciones se satistacen, 
es preferible utlizar la autocovarianza como un indí- 
cador de la textura, puesto que el promedio puede 
verse así como la componente de frecuencia cero de 
la variación espacial de los tonos de la imagen. La 
autocovarlanza representa en esta situación, la varia- 
ción de los valores de los pixels con una frécuencia 
diferente de cero. 
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La correlación cruzada de dos campos aleatorios 1 
y 9 está dada por 





lets) Er) 





onde, os la función de densidad conjunta de Hr) 
y alr): Un covarianza cruzada está definida por la 
siguiente relación 


Epia) EA) Cr lt 0,0111 


y quo so puedo roescribiren función dela correlación 
Cruzada 


Egg ye tr (e) 


Tanto la correlación como la covarianza cruzadas 
pueden ser utilizadas como modelos de similtud 
óntre dos campos aleatorios, ya que si éstos son 
gualos, los valores correspondientes de tales funcio- 
nos será máximo y su valor decrecerá a medida que 
los campos sean diferentes hasta ser completamente 
decorrelacionados, en cuyo caso 


Opglaira) 0 Do BED) + ELE IE 


Desde luego gran cuidado deberá tenerse en la 
normalización por posibles efectos en diferencias en 
geometría, escala y orientación que pudieran tener 
los campos ty g. Ahora bien, dos campos aleatorios 
ty g son ortogonales si para cualquier , yr, setiene 


Rpte) 0 


Es necesario hacer notar que lo expuesto hasta 
ahoraes válido para campos aleatorios reales puesto 
que la función de distribución se define en términos 
del número, el que se supone real. 

Un caso particular e interesante es cuando un 
¡campo aleatorio es homogéneo o estacionario; esto 
es así si se cumplen dos condiciones: a). El valor 
esperado +(1(e') . esindependiente de la posición 


r,porlo que «ir) =u = constante e independiente 
der yla segunda condición es que, b). Su función de 
autocorrelación sea invariante frente a translaciones. 


Relevo) = Ryglei 9 risas + es) 
.“ 





(592) 


onder es unvector cualquiera del plano». Como 
esta útima relación es valida para cualquier? Y os 
entonces haciendor, = 1, lafunción de autocorrela. 
ción puede escribirse de la siguiente manera. 


Rp) + Bl 950) > Kgles 02) 


5.99) 


que depende únicamente del vector diferencia, 1, 
que es finalmente una condición de estacionaldad 
espacial, En imágenes de escenas reales general. 
mente sólo puede suponerse estacionalidad local, 
aún cuando en muchas ocasiones se hacen opera: 
ciones a la imagen como si ésta cumpliso las dos 
condiciones mencionadas (5.91 y 5,92). Además de 
sto podemos decir que cuando un campo aleatorio 
+s homogéneo se tiene que 








esto quiere decir que si r, = (X,, Y;) Y fa = (%o1 Ya) 
entonces 





a y a 


ysisehace 2. 11-03-71 =1s ,la función de 
autocorrelación puede escribirse como. 





poner) + Kyg(0.8) = Rica 8) 
(5:06) 
Porotro lado como 1, es arbiario podemos esco- 
¡erto de manera conveniente, así que hacemos 1, = 
1, + 1 y entonces para dos campos aleatorios 1Y y 
cónjuntamente homogéneos se tiene lo siguiente 
Bega) = Tte e odas + 0) 
(sm 
y de aquí 


Helen e telar + 01) = Regles e rasta + 9) 


15:98) 
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¡conlo que substituyendo el valor que escogimos para 


To» tónemos. 


Bl 
palace 
Rele e) e llos — race rato) 


Ahorabien,introduciendolas definiciones para : y: 
, setiene finalmente 





que representa la condición de estacionalidad con- 
junta y que puede expresarse como una invarianza 
frenteatranslaciones. Es posible que existan campos 
aleatorios con un grado mayor de simetría, es decir 
'quela autocorrelación sea invariante ante rotaciones 
detal forma que sólo dependa de la distancia Euci- 
diana 


LLos campos aleatorios que cumplen con esta pro- 
piedad se llaman homogéneos e isotrópicos. Se 
puede afirmar en general, que cuando una imagen 
contiene estructuras más simples, el campo aleatorio 
rospectivo satisface más propiedades de simetría. 
Las condiciones dadas por (5.91 y5.92) establecenta 
estacionalidad en lo que se conoce como el sentido 
amplio, Además el concepto de campos aleatorios 
homogéneos es de hecho una generalización, a 
"espacios de dos o más dimensiones de los procesos 
estocásticos estacionarios, Otra manera de ver la 
estacionalidad espacial, es cualitativamente imagi- 
nando que la función que genera los valores de los 
pixels es la misma. Esto implica que las condiciones 
físicas bajo las cuales so adquirió la imagen son 
iguales, indepandientemente de la posición relaiva 
que guardan los CIV de la escena. Por condiciones 
físicas entendemos la geometría particular del expe- 
fimento, los requerimientos de homogeneidad de la 
fuente de radiación y la respuesta uniforme de los 
detectores del sistema de captura. Un aspecto rela- 
cionado con estos conceptos de simetría es el de la 
correlación observacionaP relativa a CIV vecinos. 
Algunos sistemas de captura funcionan detal manera 
¿que al realizar mediciones acerca de una escena lo 
hacen sobre CIV traslapados. En otros casos los 
mecanismos de degradación son tales que ciertos 
fenómenos de dispersión"! hacen que el valor de 
medida de un CIV esté infisido por susvecinos cerca- 
nos. En estos dos ejemplos se describe el problema 
dela correlación observacional cuyo principal efecto 
sel de ocasionar una falta de estacionalidad espa- 
cial en el campo aleatorio resultante. Entodo caso es 
necesario eliminar cualquier correlación observacio- 
nalatravés de métodos como la descomposición en 


componentes principales. Otra propiedad promedio 
queseempleaenla descripción de unaimagen digital 
esla que se refiero ala densidad espectral o espectro 
de potencia de un campo aleatorio homogéneo, Esta 
densidad se define como la transformada de Fourier 
dela función de autocorelación 





e | neta, hexol 2130001 + Sus Dldod8 


Obien la densidad espectral cruzada de dos campos 
aleatorios, se escribe en forma vectorial como 





Jer 


6.101) 
CONT = (1.8) yu = (usas) . Cabe anotar aquí que 
¡es muy fácil escribirla versión digital para la densidad 
espectral. Dada la intorpretación que se dio pe 
función de autocorrelación entérminos de latoxtura, 
puede ahora decirse que la densidad espectral es 
distribución de frecuencias de la variación espacial 
los valores de los pixels que componen la imagen 


digtal. 

Las propiedades promedio que se han planteado 
en esta Última sección se encuentran definidas final- 
mente atravós de la función de densidad, la que es o 
ificilo laboriosa de evaluar. Un principio que cumple 
cierta clase de campos aleatorios permite la evalua- 
«ción directa de las propiedades promedio menciona. 
das. Este principio es el que se refiero al teorema 
ergódico'” y que se verá aquí en forma cualtativa y 
relacionada únicamente con los dos primeros mo 
mentos; aún así esto es suficiente para entender el 
procedimiento de evaluación directa que se comentó 
ariba. 


En lo que sigue se supondrá que la discusión 
ú«concieme solamente a campos aleatorios homogé- 
q a 

2 Y Ralf 1), 50 que son 
relaciones paralaspropeiades promeciodeconjun 
todelcampos aleatorio, es decirson conocidas como 
promedios de conjunto. Veamos ahora qué significa 
lnexpresión paralamedia. Esta quiereclecir que 
el experimento bajo desarrollo se repite un gran 
"número de veces (los CIV medidos) y que para cada 
realización se observa la banda %., se muestrea la 
“unción bicimensional correspondiente de entre la 
famila K() en la posición y se promedian los valores 
asiobtenidos. Considérese ahora elsiguiente prome- 
do espacial 
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donde « es unaregión conectada simple acotada 
en el plano xy, 1. es elárea de tal región y por lim 
£ == se entiende un proceso que tiende a cubrir 
elplano x.y entero. Larelación (5.104) representa una 
familia de ecuaciones, una para cada banda de la 
imagen mutiespectral' $101, 

Al substituir este campo aleatorio en (5.104) se 
produce un número para cada ". . La colección de 
stos números es una función sobre el conjunto de 
bandas; de aquí que A sea una variable aleatoria en 
¡góneral. Podemos ahora decir lo siguiente: 1). Si la 
Variable aleatoria A es tal que su valor es el mismo 
paratodas las ', , es decir si A es una constante; 
Si esta constante es igual a la media — del campo 
aleatorio homogéneo, entonces éste es ergódico con 
respecto a la meda Esta conidoración de 
ergodicidad es muy. importante en aplicaciones 
computacionales relacionadas con el análisis de 
imágenes, ya que si una colección de imágenes 
(bandas) forman un campo aleatorio que sea ergódi- 
co con respecto a la media, entonces al primer 
momento de este campo puede ser obtenido por un 
simplo promedio espacial delos valores de los pixels 
decualquieradelas bandas delacolección. Consido- 
racionos análogas pueden hacerse en relación a la 
autocorrelación, la que también es un promedio de 
¡Conjunto del campo aleatorio en cuestión. La defini- 
ción para este promedio es equivalente a repetir al 
experimento respectivo un gran número de veces y 
para cada realización, observar las bandas +, 
muestrear la función bidimensional asociada a partir 
de la familia ((1) para los puntos r, y (,, determinar el 
producto de estas dos muestras y promecir el resul- 
tado así obtenido. Sea ahora el siguiente promedio 
espacial para un campo aleatorio homogéneo 








donde ; es unaregión conectada simple acotada en 
elplanoxy, ». eselárea cubierta portal región y por 
lim o +” seentiende un proceso quetiende acubrir 
elplanoxy entero, Esclaro que *(«,s) esunavariable 
aleatoria en general, ya que para cada banda se 
Obtiene un número. La colección de estos números 


para el campo aleatorio! e, 
aleatoria. Ahorabien, siesta variable aleatoria esigual 
auna constante para todo par. (2,5) y sieste valor 
constante es igual a %,,(a,s) , entonces el campo 
aleatorio es ergódico con respecto a la autocorrela- 
ción. Siun campo aleatorio es ergódico con respecto 
ala autocorrelación, entonces la función 4 (4,8) 
puede obtenerse simplemente porel promedio dado 
por(5.105) para cualquier banda del campo aleatorio. 
¡Como hemos visto.l cálculo de propiedades estacís: 
ticas de una imagen, como la media y la autocorrela- 
ción es grandemente simplicado si el proceso que 
¡generó el campo aleatorio correspondiente es ergó- 
ico. Es muy díficil demostrar que tal o cual proceso 
es ergódico con respecto a todos los momentos 
estadísticos, generalmente sólo puede hacerse para 
los primeros de ellos, Aún así, ergodicidad para los 
dos o tres primeros momentos es suficiente para 
muchas aplicaciones importantes del análisis de 
imágenes. Enlapráctica, la ergodicidad de un campo 
aleatoriose asume con frecuencia yso realizan cálcu- 
los como sital fuese el.caso. Muchas imágenes no 
satisfacen el principio de ergodicidad por lo que 
considerables errores pueden introducirse al evaluar 
Propiedades estadísticas y al realizar clasificaciones 
espectrales, Es responsabilidad del investigador el 
evaluar el error que se cometo al suponer válida la 
úergodicidad así como el de explorar su aplicación en 
'uncampo aleatorio. Podemos añadir 

úergódico es estacionario, pero un proces 
o no siempre es ergódico. También, un sistema 
ergódico es no-destructivo aunque nonecesariamen- 
to sea válido a la inversa. Un sistema con memoria 
nunca es exgódico. Para finalizar podamos decir que 
lavalidez delarelación(5.105)implicacualitativamen- 
teo siguiente: supongamos que por algún procedi. 
miento podemos obtener n-copias del sistema que 
deseamos estudiar; digamos n-muestras o n-propa- 
raciones. Dejamos ahora evolucionar en el tiempo 
estas n-copias y obtenemos imágenes en n-tiempos 
diferentes. Así por ejemplo la primera copia la obsor- 
vamos en el tiempo t,. la segunda en el t, y así 
sucesivamente. Ergodicidad quiere decir entonces 
que es lo mismo observarla evolución de un sistema 
en el tiempo que la observación de n-sistemas en 
diferentes tiempos o estados de evolución. De esta 
forma, a panir de la imagen de un bosque digamos y, 
si se cumple la ergodicidad podemos inferir el creci 
miento de un pino porla observación de varlos pinos 
sn diferentes etapas de crecimiento, 
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Introducción 


on!a intención de presentar lo más sencillamente para el lector 
el manípuleo de laimagen en una computadora, supongamos 
que la misma está dividida en una cuadrícula, donde cada 
elemento, denominado pixel, contiene un valor de gris proporcional a la 
reflectividad de la imagen en esa posición. Cuando una escena está 
compuesta por varias imágenes correspondientes a respuestas en dit 
tas bandas del espectro, su representación se hace en tres dimensiones: 
Con el fin de optimizar el 
uso de las características 
del computador, la unidad 
de almacenamiento en 
medios periféricos es por 
línea completa y por canal. 
Todos los procedimientos 
detratamiento de la imagen 
deben tener en cuenta, 
como se rescata cada 
unidad de información, lo 
queincide con gran peso en 
el tiempo de proceso total. 





























N.dela R. 


Estetrabajo pretende entregarnociones generales de procesamiento digital de imágenes y sus correspondientes ejemplos 
reales. Ha sido diseñado como una introducción al tema, especialmente útil para no especialistas. Por ell, la 
Editoialdela Rovista SELPER ha considerado de interés inciirio en este número, como una interesante iniciación altema. 
Esto trabajo fue extractado de los Apuntes Docentes del Curso 23, “Percepción Remota: Perspectivas de Avance y 
¡Cooperación Aeroespacial en Latinoamérica”, de la Temporada internacional 1986, "América Latina: Una Realidad", de la 
Universidad de Chile y dictado por los profesores Sr. Mauricio Araya F. (Universidad de Chile) y Sr. Marcelo Campi, 
de Procesamiento Digital de imágenes. 
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1. CARACTERISTICAS PARTICULARES DE LOS SISTEMAS 





ía: 


procesamiento digital de imágenes está basado en uno. varios módulos de computación. Un 
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El esquema anterior es general y por supuesto 
puede variar con cada instalación panicular. Es 
Importante destacarciertas caractorísticas especialos 
par su importancia. 


A. El uso de memoria RAM de alta velocidad 
permite alcanzar procesos en tiempo real. Como 
actualmentelatecnología digital no ofrece soluciones 
prácticas para manejar grandes volúmenes de 
información en estas memorias, se debe contar con 
medios más lentos de almacenamiento (discos) pero 
cuya capacidad de almacenamiento es muy superior. 





















































8. Elusodo un segundo computador (generalmente 
microcomputador) esclavo del computador central y 
con su propia memoria rápida (RAM), permite 
disminuirlas cargas del primero, lo que Incid 
medida en los costos finales del sistema pues se 
Puede llegar al uso de computadores personales tan 
difundidos actualmente y de un precio 
significativamente inferior a grandes 
computadoras. El principio de optimización del 
sistema es no agregar capacidades que al final nunca 
so usan pero sí inciden en el precio total, o lo que es lo 
mismo decir, se deben buscar sistemas los más 
orientados al tema. 














Ahora bien, como vimos, la memoria central es 
costosa ylimitada pero muy rápida Si éstanos permite 
almacenar varias líneas. de una imagon on la misma, 
los algoritmos deben estar condicionados a trabajar 
sobre estas líneas para producir la salida, 


De aquí podemos ver claramente que an bien de la 
velocidad del sistema, los programas del sistema 
deben ser reducidos y optimizados para realizar la 
menor cantidad de operaciones de entrada/salida a 
disco y, por el i " 
procesar, ratar de eliminar operaciones redundantes. 
durante el proceso, Un camino muy usado es usar 
tablas de conversión que son arreglos numéricos 
limitados generalmente a 256. elementos y cuya 
relación de transformación se establece unasolavez al 
principlodelprograma. Estatabla permite paraundato 
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de entrada obtener la salida pre-programada en 
forma instantánea. 


(C.El sofware del sistema debe cumplir con cienas. 
características para ser eficiente: 


3) Volocidad; Se logra optimizando algoritmos y 
manejo de datos, 





b) Confiable; Además de entregar resultados 
¿correctos debe adventr ingreso de datos erróneos u 
órdenes equivocadas. 


6) Interactivo: Debe buscar darle al usuario, 
generalmente lego en computación, la mayor 
capacidad de control a los procesos. 





d) Versátil: Permitir su crecimiento futuro. sin 
necesidad de rediseñar el sistema o romper 
"esquemas muy rígidos. 


e) Portátil: Tratar de usarlenguajes universales para 
que ante un eventual cambio de computadora central 
o sea necesario empezar desde el principio. 


2, ENSANCHES DE CONTRASTE 


La necesidad del ensanche de contraste de una 
imagen digital se basa en el principio de que la vista 
alcanza a distinguir un número finito de valores de gris 
muy inferior a los provistos por este tipo de imagen. 
Estalimitación produce confusión cuando los niveles 
de gris son parecidos. Las metodologías aplicadas 
pormiton modificar las conversiones de gris para que 
élojo humano pueda percibirtonalidades necesarias. 

Entre los tipos de ensanches de contraste 
aplicados a imágenes digitales, los más usuales son 
los lineales, y por nivelación de histograma. 


Consideremos la distribución del histograma de 
niveles de grís de una imagen: 





NO piel 


elas de 
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El ensanche de contraste lineal consiste en ele 
un límite inferior y otro superior dela curva y alítrazar 
una recta de transferencia: 








E nietos de grís 
de entrado 














De esta forma los valores de gris entre A y B se 
distribuirán entre O y 255 a la salida. 


Nivelación de Histograma: 


Consiste en asignar mayor diferenciación do grises 
en aquellos lugares más poblados de la curva 
histograma. El procedimiento es muy sencillo. 


1%) Se calcula un incremental) como el cociente 
entrela población de pixels y el número de niveles de 
¿gris deseados (p. ej. 16) 

22) Serecorre la curva contando!as poblaciones de 
grs y cada vez que se contabilizan 4, pixel, 
contabiliza un nuevo gris a la salida. 


Conversión de Datos de Entrada a Salida 
(Tablas) 











Elprocedimiento de conversión se hace através de 
tablas de conversión, también denominadas “Look- 
uptables” que permiten una rápida conversión. 
Desde el punto de vista de un programa, es 
simplemente un arreglo numérico unidimensional, 
también denominado vector, con 256 elementos 
talforma que elíndice del mismo corresponde alnivel 
de gris de la imagen original y el contenido del vector 
en esa posición es el valor de gris a ser asignado ala 
imagen transformada. 


Estetipo de procedimiento es altamente eficiente y 
¡con solo cambiar el contenido del vector se logra un 
nuevo contraste. 





Para aquellos equipos en que este procedimiento. 
de transferencia se hace a nivel hardware se puede, 
entiempo real, r moviendo un cursor hasta lograr el 
contraste deseado. 
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Aplicaciones. 


El proceso de ensanche por nivelación de 
histograma se utiliza en aquellas imágenes en que es 
necesario resaltar una subimagen bien definida 
aunque se pierda contraste en el resto de la imagen. 

Elensanche lineal pesa igualmente los tonos de la 
imagen dando un contraste más "suave" que el 
anterior usado especialmente en imágenes con 
riqueza de tonalidades. 


Pseudocolor. 


Algunos sistemas permiten colorear una imagen 
deacuerdo a unatonalidad de gris determinada. Este 
proceso, si se hace en forma interactiva (p. ej. 
marcando con un cursor en la imagen el nivel de gris 
deseado), tiene gran aplicación pues una imagen 
blanco y negro al ser coloreada, enriquece 





-Temporales 
Clasificadores de Textura. 


Son aquellos que analizan agrupamiento de pixels 
y los agrupan según características con mortologías 
matemáticas similares. Se orientan a campos que van 
desde la rugosidad de terreno hasta análisis de 
Caracteres manuscritos. 

Son además indispensables en imágenes 
obtenidas de sensores activos como el radar de 
apedura sintética, elcualcomo beneficios cuenta con 
su operatividad en todo tiempo (con nubes) y alta 
resolución en contra de su incapacidad de captar 
respuestas espectrales. 


Clasificadores Multiespectrales 
Son aplicados a imágenes que cuentan con 
respuestas muliespectrales de reflectancia. El 


principio consiste en separar los pixels dentro de un 
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enormemente la capacidad del ojo en distinguir 
tonalidades. Las aplicaciones más útiles son en 
imágenes pancromáticas (placas o imágenes 
radiográficas, por ejemplo). 


3. CLASIFICACION DE IMAGENES. 


Es un concepto del procesamiento digital de 
imágenes que agrupa los procedimientos mediante 
los cuales una computadora puede agrupar los 
elementos (píxels) de una imagen de acuerdo a 
categorías o clases. 


Normalmente se pueden agruparlosclasifcadores 
según sus características: 


-De textura 


Luego viene el procedimiento de agruparlos de 
acuerdo a un criterio determinado. Esto se hace 
mediante la selección de porciones de la imagen 
conocidas y buscando cuáles son sus ubicaciones 
dentro del espacio multiespectral, además de sus 
características estadísticas (p. ej. dispersión, matriz 
de covariancias de canales), que en conjunto forman 
la denominada "firma" espectral de la clase. 

Esto involucra un trabajo previo a la clasificación 
que se debe obtener mediante el conocimiento de la 
verdad terrestre de la clase que se desea clasificar en 
la imagen. Generalmente, en imágenes satelitarias, 
signiica trabajo de campo y a veces apoyo de 
fotografía aérea e información cartográfica. 
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Entre los métodos de discriminación más usados 
podemos nombrar: 

- Máxima verosimilitud (*Maximum llkelihood”) 
basado en la teoría de decisión de Bayes. Es uno de 
los más usados por su precisión pero a su vez es el 
más lento por la cantidad de operaciones 
involucradas, en la clasificación de cada pixelsel 

Distancia mínima: Se utiliza como criterio de 
asignación de pixellla distancia euctidea mínimaen el 
espacio multespectral entre un pixel y las distintas 
clases. 

- Hipercubos: Se basa en la inclusión de un pixel 
dentro de un hipercubo en el espacio que caracteriza 
cada clase. Produce gran cantidad de pixels con 
mezcla de clases pero es el más rápido y útil en 
aplicaciones donde no se exige gran precisión en la 
clasificación (<60%). 

Los Clasificadores multiespectrales vistos hasta 
ahora se denominan supervisados, pues existe un 
entrenamiento previo para indicar al sistema que 
clases debe buscar. 

Pero, en ciertas circunstancias, en las cuales se 
desea hacer un estudio en una zona totalmente 
desconocida, es conveniente recurrir a métodos no 
supervisados que permiten separar, en forma 
itorativa, los elementos de laimagen en clases que se 
van formando de acuerdo a crterios que pueden ser 
aportados por el usuario. El más diversificado de 
estos algoritmos es el denominado ISODATA, 

Un método no supervisado puede serunexcelente 
paso previo para encarar un trabajo de campo que 
luego desemboca en una clasificación supervisada, 




















4. MODELO DIGITAL DE TERRENO 


La incorporación del tratamiento digital de 
imágones y su uso relacionado con aspectos 
cartográficos mediante computadores, creó la 
necesidad de incorporar la información topográfica 
con formato compatible a estos métodos. 

Así, surge la creación de los modelos digitales de 
terreno, que consiste en entregar, en unidades 
geográficas. determinadas (pixels), información 
correspondiente a la altura del terreno en ese sitio. 

Dos aspectos deben ser considerados en esta 
información: 

La primera, es la forma de obtenerla por única vez 
parair formando un banco de datos geográficos yla 
Sogunda consiste en las diferentes aplicaciones que 
se pueden consegui 


4,1. Construcción del M.D.T. 


Los métodos actualmente aceptados para esta 
tarea, se basan en la incorporación de datos en una 
malla poco densay luego suextrapolaciónaunagrila 
densa a ser almacenadas digitalmente. 

La adquisición de los puntos de control, se hace 
normalmente de pares estereoscópicos en las 
escalas de trabajo o a panir de Íneas de nivel de 
mapas topográficos de la región previamente 
elaborados. El uso de mesas digitalizadoras 
incorporadas al sistema de procesamiento digital, es 


prácticamente mandatorio para obtener calidad y 
velocidad en la toma. 

Actualmente, se está trabajando en métodos de 
adquisición de cotas en forma digital y automática a 
partir de imágenes estereoscópicas digitalizadas, 
usando técnicas de correlación de patrones y 
determinando asíel desplazamiento de paralaje entre 
los pares. 


Esta tecnología, cuando pueda ser aplicada 
masivamente, ahorrará esfuerzos humanos y podrá 
ser llevada a cabo sin contar con equipos 
fotograméfricos sofisticados y de alto costo. 





42. Aplicaciones 
Estoreoscopía. 


Con el agregada de información de altura de 
terreno a una imagen se puede obtener su par 

implemente desplazando “la 
ubicación de los pixeles correspondientes 
(desplazamiento de paralaje) a cada punto de la 
imagen enuna magnitud proporcionalasu altura. Los 
resultados pueden ser observados mediante uso de 
anaglifos en un resultado color superpuesto o como 
Imágenes estereoscópicas normales. 





Imágenes Sintéticas. 


Partiendo dela basede considerar elterreno como 
una superficie lamberiana, se pueden plantear 
modelos de iluminación fjando una posición 
arbitraria del soly con el conocimiento del modelo de 
alturas, con lo cual se deduce la pendiente en cada 
pixel, se pueden generar distintas imágenes con 
sombras obtenidas del modelo de iluminación 
adoptado. 








Estaes una forma novedosa de entregar alusuario 
información topográfica sin necesidad de tomar la 
fotografía correspondiente y mejor que un plano 
topográfico que, no obstante, se puede generar lo 
mismo. 


Imágenes en perspectivas. 


Eligiendo el punto de visión, con el conocimiento 
de la alura en cada olemento de la imagen y su 
roflectancia, se puede plamear un modelo de 
proyección central que simula una vista en 
Perspectiva de la imagen venical. 








Mediante el agregado de los programas 
correspondientes se puede, por ejemplo, ejempiíñcar 
elllenado dela zona con agua (caso de unarepresa) 
y verificar visualmente qué nivel alcanzaría la represa 
sobre la imagen para distintos volámenes de agua. 


5. TRATAMIENTO RADIOMETRICO DE 
IMAGENES 


Veremos algunos delos tratamientos más usa 
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5.1. Filtrado Pasa-bajo 


El objetivo es atenuar los efectos de bordes de una 
imagen obteniendo una imagen con un aspecto 
difuso. Es útil cuando se desean disimular efectos no 
deseables en imágenes producidas por procesos 
matemáticos sobre la misma. Un caso típico es la 
continuidadde valores luego de una densificación de 
puntos por extrapolación en la construcción de un 
modelo digital de terreno. 

La forma más normal de obtener este efecto es 
generar una ir cuyos elementos son el 
promedio aritmético de los adyacentes. 


52, Filtrado Pasa-alto 


Permite resaltar los bordes obteniendo. mayor 
efecto de nitidez, El procedimiento consiste en 
amplificar el efecto enlos bordes dentro dela imagen 
(detectados por diferencias de gris entre pixels 
adyacentes que superan un umbral): 








0 


As umbral 














El eje x representa la ubicación del pixel dentro de 
una línea. 


5.3. Filtrado de gradiente 


El procedimiento, también denominado 
“pseudoplástico” permite extraer bordes de una 
Imagen. Ello se hace representando las diferencias 
entro niveles de grís como un nuevo nivel de gris. De 
acuerdo a las direcciones en que se efectúa esta 
diferencia, se obtienen diferentes efectos sobre la 
imagen de salida (ver Pág. 71 curso N* 23 
"Percepción Remota: Perspectivas de avance y 
cooperación aeroespacial en latinoamérica” y otras 
mostradas al final de este Anexo). 


5.4, Filtrado por convolución de máscaras 


Es un fito más general y que puede ser diseñado 
por el mismo operador. 





Consideremos una tabla de 9 elementos; 





co» 
zoo 
ES 











Para realizar la convolución de esta máscara, 
apoyamos la misma sobre cada cuadro de 3 líneas x 
3 pixels de la imagen normal: 


AJA Ava 
Agos Ager ajos 
Agra Merge Avajea 


| 


porel. 


——— nea). 





y realizamos la siguiente operación: 





igor Ayo Ao o 0 o Argo do 
Aajor 

Donde B,, es la ubicación del pixel en la línea 1 y 
pixel k de la imagen de salida. 


De esta forma el usuario puedo probar distintos 
coelicientes [ a, b...... 1]y Obtener con cada uno de 
ellos un trado ciferento. 

Esta metodología es especial para optimizar anivel 
del computador de representación visual pues, al 
tener éste acceso a la memoria RAM permito hacer 
"operación en un tiempo equivalente alde refresco de 
Pantalla (denominado para estos casos tiempo real). 

SÍ, además se generan los coeficientes (a,b....ila. 
partir de una coordenadalx,y del cursor, se obtiene 
un efecto mediante el cual el usuario moviendo el 
cursor va viendo los distintos ftrados de la imagen 
hasta seleccionar la configuración Óptima que luego 
esaplicadaa la imagen completa, de mayorvolumen, 
que reside en disco. 











5.5. Transformación de componentes 


principal 


Supongamos que representamos los pixeles de la 
imagen en un espacio multiespectral. 

Si consideramos que las distribuciones de gris en 
cada canal son gaussianas, su representación en un 
espacio mullespectral será un elipsoide. 

El procedimiento es crear un nuevo sistema de 
coordenadas a parir de una transformación del 
original detal forma que el canal 1 del nuevo sistema 
coincida con el eje mayor del elipsoide. 
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De esta manera, la proyección delos pixeles sobre 
este canal, por ser el ejo mayor de la distribución, 
presentará la mayor dispersión de valores de gris, 
dando una mayor discriminación de características 
dela imagen. 

Por supuesto, las imágenes obtenidas sobre estas 
proyeccionesnorepresentanvalores dereñectividad, 
'más bien contienen una integración Óptima de los 
valores de gris encadacanaldetal manera de obtener 
Untesumen de la imagen multiespectral ena primera 
componente principal, degradándose hasta teneren 
la úlima componente valores de gris no 
úcorrelacionados. 





















6. Complemento 


A continuación se presentan varios ejemplos 
fotográficos que, aunque están en blanco/negro, dan 








Imagen Landsat sobre una zona de Santiago de Chile, procesada computacionalmente en el Centro de 
Sensores Remotos dela Comisión Nacional de Investigaciones Espaciales (CNIE) deArgentina, en programas 
de cooperación con la Universidad de Chile. La imagen de arriba corresponde a un proceso especial, en que 
se ha determinado un modelo de cotas para producir una imagen abatida o en perspectiva. 
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E 1 + PROCESADORA LANDSAT CNIE - ARGENTINA + 
1 2-247-076 4 B1 1143-1800 A 11DEC82 


Ejemplos de productos procesados computacionalmente en el Sistema Pl del Centro de Sensores Remotos 
de la CNIE, Argentina 
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PRINCIPALES ACTIVIDADES 
DE SELPER (1980-89) 


PRIMERA REUNION PLENARIA SELPER. Quo, Ecuador no- 
'embre24-27, 1980. Cone princpalauspicio nacional gel Insisto, 
'Goográo Nacional y CURSEN y el apoyo ntemaciona! de U.S. 
AQUNOy lr Intomational Development 


SEGUNDA REUNION PLENARIA SELPER Y TALLER DE TEC- 
MOLOGIA DEL SATELITE SPOT. Quio, Ecuador, abel 1216, 
1982. Con el principal auspicio naciona de ION y CURSEN y el 

ntemacional del Centre Nañonas O'Ensdes Spatales 
(CES) on Francia. 


TERCERA REUMION PLENARIA SELPER Y TALLER CANA- 
DIENSEPERCEPCION REMOTA. Santago, Che, noviemire21- 
25, 1983 (Auditorium IDIEM, Facutad de Ciencias Fisicas y Malo 
máticas, Universidad de Che), Con el auspicio. nacional dela 
Universidad de Chile (Departamento de Investigación y Bmtotecas 
y Depto. Geología y Geofísica) y el pairocino interacional de 
“Canada inoraional Development Agency (CIDA). 


CUARTA REUNION PLENARIA SELPER. Santago, Chio, no. 
Vembre 1236, 1984. (Aucitonum IDIEM, Facunaa de Ciencias 
Físicas y Matemáticas, Unwersidad de Cn) Con el peneipal 
auspicio, Intemacional del Centre National D'Eludes Spatales 
(CNES) de Francia, Canadá Centre for Remoto Sonsing (CCA) 
Embajadas de España y Financia, ONU/FAO, ONU/CEPAL 
Commitos on space Rasearch/interational Counea ot Sientc 
union (COSPAR/ICSU), ntemational Union of Geología: Scien- 
cios (UG) y otras agencias. 


QUINTA REUNION PLENARIA SELPER Y SIMPOSIO CANADA/ 
IBEROAMERICA DE PERCEPCION REMOTA, Otiana, Canadá, 
julio 15-19, 1985. Con el principal auspicio dl Miistao de sta: 
lonas Exteriores de Canadá y a colacoración del Centro Cana- 
dlense de Porcepción Remota (CAS); comaó, organizador 
'SELPER (Sede Contra en Unwersidad de Chio, Depto. Geología y 
Geofísica. 


TALLER DE GEOLOGIA REGIONAL UNESCO/1UGS/SELPER. 
¡Seo Josh dos Campos INPE, Brasi, diembre 2-4, 1985, Con el 
principal auspicio de UNESCOe ntémañional Union l Geologia! 
¡Sciences (UGS), para coorainar programas. Geológicos Regiona- 
os on Latinoaménica (Zona Amazónica y Andina). con SELPER 
como organización aglutinante y coordinadora. 


CURSO DE PERCEPCION REMOTA Y PERSPECTIVAS DE 
AVANCE Y COOPERACION AEROESPACIAL EN LATINOAME- 


SIMPOSIO NACIONAL SOBRE LA EVALUACION DE LA ERA 
ESPACIAL EMMEXICO, México D.F. México, 15 19e mayo 1988, 
"Organizado por GIAE/UNAM. SELPER,Inutaco Oti 


REUNION DE NACIONES UNIDAS SOBRE COOPERACION 
REGIONAL EN CIENCIA Y TECNOLOGIA ESPACIAL Y SUS 
APLICACIONES EN AMERICA LATINA. Cartagena de Indias. 
Colombia, Junio 162, 1988. Organizado por ONU (Didsión ds 
Espacio Exteior_y CEPAL) y patrocinada por el Gotiemo de 
Colombia. SELPER, Inutado Oficial 


00Vl REUNION PLENARIA COSPAR/ICSY (Commitiee on Spa: 
ca esearch/internallonal Counct ef Sesent: Unions) Toulouse, 
Frencia,juvo 30 -Julo 12, 1986, SELPER, Co patrocinante del 
"TalorX Percepción Remota de ntorés para Prísesen Desarroso”. 
También co patrocinado por ONU/FAO, ONU/COSTED, WMO, 
WHO, COSPAR/ICSU, EARSEL y otras agencias. 


SEXTA REUMION PLENARIA SELPER - SIMPOSIO LATINOA- 
MERICANO DE PERCEPCION REMOTA -1V SIMPOSIO BRASI- 
'LENO DE PERCEPCIÓN REMOTA. Gramado, RS, rasa, agosto, 
10.5, 1988. Reunión realizada en cooperación con el nsiuto e 
PesquisasEspacils (NPE) la Sociedade Brasisrade Cartograa. 
Geodesia, Fotogrametría e Sensoriamento Remoto (SEC) y otras 
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MAIN SELPER ACTIVITIES 
(1980-89) 


FIRST PLENARY MEETING OF SELPER,Qulo, Ecuador, novem: 
úDet24:27, 1980, Under ne national suppor ol CLIASEN and IGN 
ana the intemañonal Suppor 0l U.S. Agen for Intemational 
Development (USADO) 


SECOND PLENARY MEETING OF SELPER AND WORKSMOP 
ON SPOT TECMNOLOGY. Quito, Ecuador, apri 12:10, 1982, 
Uncer ne national support CLIRSEN ana IGN. and ne ntoma: 
úonalsuppor ol Centre National D'Etudes Spañalos (CNES) ano 
SOTA from France 


THIRD PLEMARY MEETINO OF SELPER AND WORKSHOP ON 
CANADIAN REMOTE SENSING TECHNOLOGY Santiago. Crue 
orember21-25, 1983. Held lo Universidad de Che, under Ino 
main national suppor of Untwersidad de Chio (DIB, Depa 
Geología y Gootsica) SAS. INACH and the international suppos 
Canada Intemalional Development Agency (CIDA). 


FOURTM PLENARY MEETING SELPER, Sanilago, Ch, novern 
bar 12-18, 1984. Hela sto Untorsidad de Cil, under ino main 
úañonal suppon cf Unwwersidad de Chile (DIB, Geología y 
¡Geofíca), SAF, INACH ana Iho intornainasuppor ot several 
"agencies suchas Centre Nadia D'Eludes Spailes (ONES) tom 
France, Canada Embassy, Spain Embassy, UN/ECLA. UN/FAO, 
Comites on Space Research/Iiemalional Counel ol Sent 
Unicas (COSPAR/ICSU),Intemational Union cl Geological Sc 
únces (UGS), Insituto de Posquisas Espaciis (NPE) from Brazil 
an ers. 


FIETH PLENARY MEETING OF SELPER ANO CANADA/LATI- 
NAMERICAN SYMPOSIUM ON REMOTE SENSING. Oti, 
Canada, jUy 15-19, 1985, Under 1o main sponsoring and organ 
zatonciDepsrtmert ol Eternal Atar o Canada ancho coopera. 
on ol Canada Centre lor Remote Sensing (CCRS). Important 
sessions ol Ino 1UGS ana COSPAR/ICSU were also developed al 
SELPER sessions. 


"WORKSHOP ON REGIONAL GEOLOGY AND REMOTE SENS- 
NG. Sao Joséidos Campos, Braz december 2-4, 1985, Undoriha 
main ntemmañonal sponsoring ol 1UGS and UNESCO, Gars Pro 
¡ram (Geologia! Appticaons ol omoto Sensing).sponsoreg by 
IUGS/UNESCO. 


COURSE ON REMOTE SENSING AND ADVANCES AND PER- 
SPECTIVES OF COOPERATION IN LATINAMERICA. interna: 
Bonal Season ol Unwersay of Che (Latin America: A Realty). Co. 
ponsorea by SELPER Imowgh professor ana importan Docent 
Documents. Santiago, Cri, January 13-24,1988 


NATIONAL SYMPOSIUMABOUT EVALUATION OFYHE SPACE 
AGEINMEXICO. MEXICO,D.F. México, may 15-19, 1986. Organ. 
(200 by: GIAE/UNAM. SELPER. oficia inte. 


UNITED NATIONS MEETING ON COOPERATION AND APPLI- 
CATIOM OF SPACE SCIENCE AND TECHNOLOGY IN LAT 
NAMERICA. Cartagena de incias, Colombia, Juno 18-22 1986, 
Organized ty Unted Nations (Quler Space Atar ad ECLA) and. 
sponsored by Ine Colombian Govemmoni. SELPER afcaly in- 
vea. 


LO0VIPLENARY MEETING OF COSPAR/ICSU Toulouse, France, 
June SOjuy 12, 1986, SELPER cosponsoreg Ino Worksnop X 
REMOTE SENSING OF INTEREST FOR DEVELOPING COUN: 
TRES”. cosponsorea also by UN/FAO, UN/COSTED, WMO, 
"WO, COSPAR/ICSU, EARSEL nd oners agencies 


SIXTM PLEMARY MEETING OF SELPER-LATINAMERICAN 
SYMPOSIUM ON REMOTE SENSINGAV BRAZILIAN SYMPO- 
SIMON REMOTE SENSING. Gramado, RS, Braz august 10-15, 
1986. Meeting camed Cu ná Ine cooperation ol Insituto do. 


A e o, Sep cie 1986 SELPER 


“agencias nacionales  intomacionalos 


CONFERENCIA INTERMACIONAL SOBRE NECESIDADES DE. 
ENTRENAMIENTO Y ASISTENCIA A LOS USUAMIOS EN PAL- 
SES EN DESARROLLO. Bern Occidental sepllembro 1-0, 1986, 
Founión promovida porel ntemational Economic Summit, organi 
zada pora Funcación Alemana para el Desarrolo Intemacional 
(OSE) on la cual SELPER participó en calidad de entidad invtada. 
cialmonto 


REUNIONDELAS NACIONES UNIDAS (ONU) SOBRE CiENCIA 
Y TECNOLOGIA ESPACIALES Y SUS APLICACIONES DENTRO. 
DEL MARCO DE LOS SISTEMAS DE ENSEÑANZA, Co Palio 
ada por olGobiemo de Móxico. México, D.F, Máxico, Ocubro 13. 
17, 1968, SELPE invilado Onil, 


PRIMERA REUNION LATINOAMERICANA SODRE SATELITES. 
DE AUSQUEDA Y RESCATE Sao Jonb dos Campos, SP, Bra. 
úctubre29-30, 1989, Organizada por INPE y apoyacia por SELPER. 
ara ol Programa SARSAT (Saló1os do Búsqueda y Roscalo) 
esarolado por Canada, EE.UU, Francia y UNS, 


SEGUNDO SEMINARIO NACIONAL DE PERCEPCION REMO: 
TAY SUS APLICACIONES. Organizado porla Asociación peruana, 
o Abrololografa Aplicada (APAFA) Co Palfocinada por SELPER. 
Lima, Poró, Mayo 11:15, 1997, 


CURSO DE ENTRENAMIENTO SOBRE EL SATELITE DE RA- 
DAR EUROPEO (ERG-1). Organizado y Co patrocinado por la 
'Agoocia Espacial Europea (ESA) con SELPER. Buenos Aros, 
gonna: Novtombro 4:12, 1087, 





SEPTIMA REUNION PLENARIA SELPER, SEGUNDO SIMPO- 
'SIO LATINOAMERICANO DE PERCEPCIÓN REMOTA, ll IA 
POSIO COLOMMANO DE PERCEPCION REMOTA, Coloquio 
Intornacional SPOT, Primera Rounión JOBRESA (Joint aura ol 
Fomolo Sansing Actis). Organizado por SELPCR, nllulo 
Geográtco Nacional Agustin Codazz1 (IGAC) y Sociedad Colom. 
lana de Fotogramolra y Porcapeión Remota (SC), con el apoyo. 
de avarsas organtzaciones intemacionalos (ESA, CNES, SPOT 
IMAGE, DORNIER. INPE y oras), Bogol, Colombia. Novembre 
1620,4967, 





CURSO DE ENTRENAMIENTO SOBRE EL SATELITE DE RA: 
DAR KUROPRO (ERS-1). Organizado por la Agencia Espacll 
Europos (ESA) como colaboración a SELPER. Bola Calforia. 
México. Noviombro 1900 


OCTAVA REUMION PLENARIA SELPER Y TERCER SIMPOSIO 
LATINOAMERICANO DE SENSORIAMENTO REMOTO. Orgari 
zada por SELPER yla Unworsidad Nacional Autónoma de Ménco. 
(UNAM) con elapoyo de diversas organizaciones internacionales 
(on ospacial EOSAT), Acapulco, MÓXIco,dciembro 59, 1080, 


REUNION DE EXPERTOS PARA REVISION FINAL DEL DICCIO: 
NARIO SELPER DE PERCEPCIÓN REMOTA Patrocinada por ia 
Agencia Canacionsa para Desarrol intoracional (ACON) que 
“donará SELPER este Diccionari, producto del Proyecta PERCEP 
CaaóPorú. Pachacamas (Lima), Perú, Junio 5:10, 1089 


NOVENA REUMION PLENARIA SELPER, CUARTO SIMPOSIO 
LATINOAMERICANO DE PERCEPCION REMOTA, CURSO 
INTERNACIONAL SPOT Y OTRAS ACTIVIDADES ASOCIADAS. 
Organizada por SELPER. Insituto de Tecnología Industrial (TI, 
Comisión Nacional do Investgacionos Espacilos (CNE), con el 
apoyo de diversas organizaciones Intormaciónalos. Barvoc 
Agontina, noviembre 2024, 1989. 


OTRAS ACTIVIDADES DIVERSAS PROMOVIDAS POR LA 
SEDE SELPER/BRASIL, Incluyondo co-organización y co patro 
nio e Cursos, Seminarios, Congresos, quo serán defiadas an 
las primeras publicaciones de la Hovisia SELPE en 1990. 








Pesquisas  Espacials: (NPE), Ino Sociedade: Brasileira de 
Cartograf, Geodesia, Fologramotría e Sensoriamonto Remoto. 
($80) ana other national and ntomational agencies, 





INTERNATIONAL MEETINO OM TRAININQ NEEDS AND USER 
ASSISTANCE. Barín, AFA, augus!31-soplomber 6, 1986, Organ: 
izo by the German Foundation for ntomational Dovolopmont 
(0SE). SELPER parbcipaton as Invited socay. 


UNITED NATIONS MEETING ON SPACE SCIENCE AND TECH: 
NOLOQY IN THE EDUCATION FRAME. Co-spontored by Max 
can Gowemment. Mexdeo D.F, Moxco. October 19-17, 1985, 
SELPEA oficia ina. 


FIRSTLATINAMERICAN MEETING ON SEARCH AND RESCUE 
SATELUTES. Sao Josú dos Campos, SP, Aral Octonor 20:30, 
1996, Organizad by INPE (Graslan Spaco heseareh nstuto) and 
cosponsorod by SELPER. lor no SARSAT Systom developed by 
Canada, Franco, USA, USSA. 


SECOND NATIONAL SEMINAN ON REMOTE SENSIMA AND 
TS APPLICATIONS, Organizod by Peruvian Soctty ol Prato 
¡ramenatoy and, homote Sensing (APAFA) and co-sponsored by 
SELPER. Lima, Por May. 11-15 1987, 


TRAINING CQUASE OM EUROPEAN RADAR SATELLITE (ENS+ 
1), Organized y European Spaco Agoney (ESA) andco-sponsored 
by SELPER Buenos Atos, Argentina. NoVombor 4-12, 1997, 


SEVENTH PLENARY MEETING OF SELPER, SECOND LATI- 
NAMERICAN SYMPOSIUM ON REMOTE SENSINA, THIRD CO; 
LOMBIAN SYMPOSIUM ON REMOTE SENSING, INTERNA: 
TIONAL SPOT COLLOQUIUM, FIRST MEETINO OF JOBNESA 
(ont Boara ol Fomoto Sonsing Actos), Co-organizod by 
'SELPER Agustin Codazzl National Gaographicinstute (GAC) aná 
Colombiam Soc. ol Protogrammetry and Remoto Sensing 
(SC), Cosponsored by several Intoralional organizanons. 
Bogota, Colombia. Novembor 16-20, 1907, 





TRAINING COURSE ON EUROPEAN RADAR SATELLITE (ENS: 
1), Organizod by European Space Agency (ESA) lo Indorso 
SELPER actas, Baja Calfomia, Moxico. November 1908. 


IOMTM PLENARY MEETING OF SELPER ANO THIRD LATI- 
AMERICAN SYMPOSIUM ON REMOTE SENSINA, Cosorgan. 
1209 by SELPER and National Univoraty of Mexico (UNAM). Co: 
úponsorOg by soveal Intomalional organizations. ACapuico, Mex. 
leo, Decombor54, 1988. 





MEETING OF LATINAMERICAN EXPERTISES FOR FINAL RE- 
VIEWINQ OF THE REMOTE SENSIMO DICTIONARY OF 
SELPER. Sponsored by Canada International, Development 
Agony (CIDA) and cosorganizod by SELPER, ONERN. COMIDA. 
Pachacamac (Lima), Poru Juno 5:16, 1980, 


NINTM PLENARY MEETING OF SELPER, FOURTH LATI- 
NAMERICAN SIMPOSIUM ON REMOTE SENSING. AND ASSO- 
CIATED ACTIVITIES, INCLUDING INTERNATIONAL COURSE 
ON SPOT TECHNOLOGY, Co-organizod by SELPER, INTI (na- 
"ona! Iso or Industrial Tectinology), UNLU (Lujan University), 
NIE (Argontinian Commission for Space Resoarch), Cospon: 
010 by several Internallonal organizations, Bartocho, Argentina. 
November 2024, 1989. 


MORE IMPORTANT ACTIVITIES PROMOTED BY SELPER 
SEAT BRAZIL, including co organizalon and co sponsoring ol 
Courses, Seminars, Congress insido ana ulsidoBrazl wc wal 
bo Inclucod n next SELPER Moviwes (190). 
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APOYO IBEROAMERICANO / LATINAMERICAN SPONSORING 





SELPER agradece atodos los organismos e instituciones iberoamericanas 
quehan apoyado permanentemente sus actividades, yasea sirviendocomo 
sede a sus Capítulos Nacionales o en la organización de simposios y otros 
eventos u otras actividades desarrolladas por SELPER en estos años. 
Como reconocimiento, se incluyen aquí las principales entidades 
colaboradoras, esperando contar con su valiosa cooperación en el futuro. 





ARGENTINA 

COMISION NACIONAL DE INVESTIGACIONES 

ESPACIALES (CNEL 

JTITUTO DE TECHOLOOIA INDUSTRIAL. CINTI) 
OL INVESTIGACIONES APLICADAS 





us 
DE CIENCIAS ESPACIALES (LIACE 
AEROLINEAS. ARGENTINAS 

UNIVERSIDAD DE LUJAN, UNIVERSIDAD DE 
LA RIOJA, SVDSA, AACE 





erase 
INSTITUTO DE PESQUISAS ESPACIANS 

a) 

SOGIFDADE: BRASILEIRA DE FOTOGRAMETRIA E 
SENSORIAMENTO REMOTO (5AC) 

AEROSII, — COOEVASF, — PROSPEC, HOTEL 
SERRANO, DP MINERACAO, SENSORA, EMIRAER, 
COMOMA, ANEA, AVIBRAS, 180F, FAFESP, 
AAEROFOTO CRUZEIRO,  TERRAFOTO, VARIO 


ome 
UNIVERSIDAD DE CHILE 

INSTITUTO ANTARFICO CHILENO: (1NACI) 
INSTITUTO GEOGRAF ICO. MILITAR. 10M) 
SERVICIO. AEROFOTOGRAMETR ICO. (SAF /EACW) 
EXPRESA CONSTRUCTORA DEL. CHOAPA LIDA. 
EOITORIAL UNIVERSITARIA 

GODELCO, SAG, INFOR, LAN CMILE 


coLoMBIA 
INSTITUTO GEOGKAFICO- AGUSTIN CODAZZI 

IGAC y oxciAr) 

SOCIEDAD: COLOMBIANA DE FOTOGRAMENRIA Y 
PERCEPCION REMOTA. (SCF) 

ASOCIACIÓN DE CAFETEROS DE COLOMBIA 
AVIANCA 





ECUADOR: 
CENTRO DE LEVANTAMIENTOS INTEGRADOS DE 
RECURSOS NATURALES POR SENSORES REMOTOS. 
COLISEO 





ECUADOR (cont) 
INSTITUTO. GEOGRAF ICO. MILITAR. DEL. ECUADOR 
ECUATORIANA DE AVIACION 





mexico, 
UNIVERSIDAD NACIONAL AUTONOMA DE MEXICO. 
CUNA 

INSTITUTO DE GEOGRAFIA: Y ESTADISTICA 
neo 


pen 
OFICINA NACIONAL DE EVALUACION DE 
RECURSOS. NATURALES. (ONE 

COMISIÓN NACIONAL. DE. INVESTIGACION: Y 
DESARROLLO AEROESPACIAL. (CONIOA) 
ASOCIACION PERUANA DE. AEROFOTOGRAF 14 
APLICADA (APAFA) 

AEROPEAU 

INSTITUTO GEoF1S1CO VEL. PERU 


CAPITULOS. SELPER EN FORMACION 


BOLIVIA 
SERVICIO GEOLOOICO BOLIVIANO: (GCODOL) 
COSTA RICA 

INSTITUTO, GEOGRAF ICO NACIONAL £10N) 
EL SALVADOR 

MINISTERIO DE AGRICULTURA Y GANADERIA. 
GUATEMALA. 

INSTITUTO: GEOGRAFICO MILITAR 

PANA. 


UNIVERSIDAD: NACIONAL: DE PANAMA. 
REPUBLICA. DOMINICANA 

INSTITUTO GEOGRAFICO UNIVERSITARIO 
VENEZUELA 

:OCIEDAD VENEZOLANA DE FOTOGRAMEIRIA Y 
PERCEPCION REMOTA. 

PARAGUAY 

INSTITUTO GEOORAF ICO: NACIONAL. 

URUCUAY 

MINISTERIO DE AGRICULTURA. Y PESGA 





SELPER thanks to all the Iberoamerican Agencies and Institutions that 
supported SELPER permanently either being the branch of the National 
Chapters or in the organization of Symposia or other events and activities 
developed by SELPER. In recognition, the main supporting entitites are 
being included below or in additional pages. We also hope your valuable 
cooperation in the future. 
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APOYO INTERNACIONAL / INTERNATIONAL SPONSORING 





SELPER agradece a todas las empresas y organizaciones internacionales 
que han apoyado sus actividades en el pasado, esperando poder contar 


ELPER 


con su valiosa cooperación en el futuro. Como reconocimiento, se incluyen 
las principales entidades colaborantes y, en páginas especiales, aquellas 
jue más se han destacado en el año correspondiente a esta Revista 





ALEMANIA FEDERAL. 
AGENCIA ALEMANA PARA EL. DESARROLLO. 
INTERNACIONAL, (DSE) 

INSTITUTO DE. INVESTIGACIONES 
AEROESPACIALES (OF VAR) 

DORNIER, ZEISS, MBB, LUFTHANSA 


CANADA. 
AGENCIA CANADIENSE PARA EL. DESARROLLO 
INTERNACIONAL CCIDA) 

MINISTERIO DE RELACIONES EXTERIORES 
DEL. CANADA 

CENTRO CANADIENSE DE PERCEPCION REMOTA 
(CCRS) 

INTERA, MONITEO, DIGIM, MDA, PAMAP 

CP AIR 


ESPAÑA 

MINISTERIO DE TRANSPORTES, TURISMO Y 
COMUNICACIONES 

EMEBAJADA DE ESPAÑA EN CHILE 


ESTADOS UNIDOS 
UNIVERSIDAD DE MURRAY (MARC) 

NASA, USGS, ERIM 

EARTH OBSERVATION SATELLITE COMPANY 
KEOSAT) 


FINLANDIA 


VAISALA (Moteorologfa) 


FRANCIA 
CENTRE NATIONAL D*ETUDES SPATIALES 
(ONES) 


FRANCIA (cont») 

SOCIETE EUROPEENE DE PROPULSIÓN (SEP) 
SPOT IMAGE 

G.DoTeA» 

AIR FRANCE 


HOLANDA 
INTERNATIONAL. INSTITUTE FOR: AEROSPACE 
SURVEY AND EARTH SCIENCES (ITC) 


PORTUGAL. 


AGENCIAS. INTERNACIONALES 


EUROPEAN SPACE AGENCY (ESA) 


NACIONES UNIDAS (DIVISION DEL ESPACIO: 
EXTERIOR, FAO, UNEP, COSTED, CEPAL, 
UNESCO) 


COMMITTEE ON SPACE RESEARCH (COSPAR)/ 
INTERNATIONAL. COUNCIL OF. SCIENTIFIC 
LNIOS (1CSU) 


INTERNATIONAL. UNION: OF. GEOLOGICAL 
SCIENCES (1UGS) 


INTERNATIONAL SOCIETY OF PHOTOGRAMMETAY 
ANO. REMOTE SENSING: L1SPRS) 


AMERICAN SOCIETY OF. PHOTOGRAMMETRY- ANO 
REMOTE SENSINO (ASPRS) 


INTERNATIONAL. ASTRONAUTICAL. FEDERATION 
0) 





SELPER thanksto all the International Companies and Organizations forthe 
valuable support givento the activities of SELPER in the pastand hopes that 
will also cooperate in the future. In recognition, the main supporting entities 
are being included below. In special additional pages, we are including the 
most outstanding entities that supported SELPER activities specially during 
the year corresponding to this SELPER Review. 











PER sepentre 1908 


Earth Observatlon Satelllte Company 





Septiembre 1988 SELPER 
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La misión fundamental del Instituto. 
Geográfico Militar es efectuar la 
representación gráfica del territorio nacional 
en mapas, carias geográficas y planos 
topográficos, 

Con tal objeto, el Instituto realiza trabajos 
de terreno y de gabinete; los primeros 
están destinados a obtener valores de 
coordenadas y elevaciones de puntos, 
mediante mediciones astronómicas, 
geodésicas, topográficas, gravimétricas y 
'geomagnélicas. Los segundos están 
destinados a la planificación de 

proyectos, al cálculo y compensación de 
coordenadas, a la restitución fotogramétrica, 
al dibujo cartográfico y a la impresión 
offset 

Los resultados de sus mediciones se 
utlizan también en los estudios de las 
ciencias de la Tierra y especialmente en la 
investigación que hace la Geodesia sobre 
la forma y dimensiones de nuestro planeta. 






INSTITUTO GEOGRAFICO MILITAR 
Nueva Santa Isabel 1640 — Santiago-Chile 


INSTITUTO GEOGRÁFICO MILITAR 
Nueva Sta, Isabel N' 1640 Fono: 6968221-8 














PRESIDENTE SELPER: 
ES 

AN 
OP. 515, Sao José dos Campos, SP, Brasil 
F% (0123) 229977, TLX (0123)) a530 INPEBR: 
'SEDE ECUADOR 


SEDE CHILE 
1980-1983 1883-1988 





DIRECTOR EDITORIAL SELPER: 

PROF. MAURICIO ARAYA FIGUEROA. 
Presidente Madero 795 (Nuñoa). Santiago, Chile 
E: (0562) 6967883, TLX (0562) 240501 BOOTH CL. 
FAX (0562) 6981474, 


¡SEDE ARGENTINA SEDE PERU 
1989-1991 1991-19891 


IMPRESO EN LOS TALLERES DEL INSTITUTO GEOGRAFICO MILITAR DE CHILE 
Nueva Santa isabel 1640. Santiago, Chile. F: (0562) 6968221, TUX: (0582) 441677. FAX: (0562) 8988278, 








